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SECTION  I 
INTRODUCTION 


A. 


THE  SEISMIC  RESEARCH  OBSERVATORY  SYSTEM 


It  has  been  noted  (Robinson,  1967)  that  seismic  data  recorded 
by  surface-sited  instruments  are  degraded  or  obscured  to  a significant  degree 
by  wind-induced  earth  tilts.  Theoretical  data  and  tests  indicated  that  this 
wind-induced  noise  component  decreases  rapidly  with  depth.  Therefore,  it 
is  possible  to  enhance  the  seismic  data  quality  by  locating  the  seismic  sensor 
at  a depth  of  approximately  100  meters.  The  Seismic  Research  Observatories 
form  the  first  network  of  single  sites  to  implement  this  observation. 


The  Seismic  Research  Observatories  data  acquisition  and  re- 
cording system  has  been  described  in  detail  previously  by  Strauss  (1976). 
Therefore,  only  a brief  description  of  the  system  is  presented  here. 


Broadband  seismic  energy  is  recorded  by  force-balance  type 
seismometers  which  produce  an  output  proportional  to  earth  acceleration  over 
the  frequency  range  of  0.02  to  1.0  Hz.  Both  long-period  and  short-period  data 
are  produced  from  each  sensor  by  selectively  filtering  the  broadband  output. 
The  long-period  data  are  digitized  and  recorded  continuously  on  an  800  bit  per 
inch  magnetic  tape.  The  vertical  component  of  short-period  motion  is  input 
to  an  automatic  detector,  which  permits  recording  of  vertical-component 
short-period  motion  only  when  certain  operator-specified  conditions  are  sat- 
isfied. Each  magnetic  tape  can  hold  up  to  fourteen  days  of  digital  seismic 
d.  a.  When  full,  these  magnetic  tapes  are  shipped  to  the  Albuquerque  Seis- 
mological  Laboratory  at  Kirtland  Air  Force  Base  in  New  Mexico.  Copies  of 
these  tapes  are  sent  to  the  Seismic  Data  Analysis  Center  at  Alexandria, 
Virginia. 
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Long -period  data  are  sampled  once  per  second.  The  instrument 
response  peaks  at  a period  of  25  seconds  with  a quantization  factor  of  5 com- 
puter counts  per  millimicron  of  ground  motion. 

Short-period  data  are  sampled  20  times  per  second.  The  instru- 
ment response  peaks  at  a period  of  1 second.  The  1 second  quantization 
factor  is  2000  computer  counts  per  millimicron  of  ground  motion  with  the 
following  exceptions.  Beginning  1 May  1976  at  Guam,  14  April  1976  at 
Wellington,  New  Zealand,  and  13  May  1976  at  Taipei,  Taiwan,  short-period 
data  were  quantized  at  2 computer  counts  per  millimicron  of  ground  motion 
to  prevent  data  clipping. 

Also  evaluated  in  this  report  are  the  Abbreviated  Seismic  Re- 
search Observatories  (ASRO).  These  stations  feature  surface-vault  seismo- 
meters rather  than  the  borehole  type.  Data  sampling  rates  are  identical  and 
response  characteristics  are  very  similar  to  those  of  the  SRO  sites.  Nor- 
malized response  characteristics  for  ASRO  and  SRO  instruments  are  shown 
in  Figure  1-1.  Quantization  factors  at  1 and  25  second  periods  are  10  and 
1000  computer  counts  per  millimicron,  respectively. 

B.  THE  EVALUATION  TASK 

This  report  presents  the  results  of  a continued  evaluation  of 
the  Seismic  Research  Observatories.  The  specific  goals  of  this  evaluation  are: 

• To  estimate  the  data  quality  and  reliability  of  each  station 

• To  investigate  the  short-period  and  long-period  noise  field 
characteristics  of  each  station 

• To  evaluate  the  performance  of  the  short-period  automatic 
detector 

• To  estimate  the  detection  capability  of  each  station 
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NORMALIZED  INSTRUMENT  RESPONSE  CURVES 


To  estimate  the  discrimination  capability  of  each  station  on  a 
regionalized  basis. 


These  evaluation  goals  are  presented  in  the  following  manner. 
First,  a suitable  suite  of  seismic  events  and  noise  samples  was  constructed 
from  available  event  lists.  The  procedure  for  selecting  this  event  suite  is 
described  in  Section  II.  This  section  also  describes  the  manner  in  which  the 
events  and  noise  samples  were  processed  and  gives  estimates  derived  from 
processing  experience  of  the  quality  of  the  Seismic  Research  Observatory 
data  and  reliability  of  the  individual  stations.  Section  III  discusses  the  short- 
period  automatic  detector  resident  at  each  station.  The  discussion  of  this 
topic  covers  the  present  design  and  evaluation  of  the  automatic  detector. 

In  Section  IV,  the  local  noise  field  characteristics  at  each 
evaluated  SRO  station  are  presented.  The  local  noise  field  is  characterized 
by  the  RMS  noise  levels,  monthly  RMS  noise  level  trends,  p^ak  noise  ampli- 
tudes, and  spectral  content  of  the  noise.  Section  V presents  the  detection 
capability  of  each  station  in  terms  of  the  entire  area  of  interest  and  in  terms 
of  specified  regions  within  the  area  of  interest.  In  the  course  of  estimating 
detection  capability,  the  effects  of  mixed  events  and  system  malfunctions 
on  these  estimates  are  also  discussed.  Section  VI  considers  the  question  of 
discriminating  between  earthquakes  and  presumed  nuclear  explosions.  The 
technique  used  in  this  section  is  the  M^-  m^  discriminant.  Section  VII  sum- 
marizes the  results  of  this  evaluation  and  suggests  future  avenues  of  inves- 
tigation to  improve  these  results.  Section  VIII  lists  the  references  cited  in 
the  text  of  the  report. 
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SECTION  II 
THE  DATA  BASE 


A.  DATA  AVAILABILITY 

The  Seismic  Research  Observatory  network,  when  fully  opera- 
tional, will  consist  of  at  least  the  eighteen  stations  listed  in  Table  II- 1 and 
shown  in  Figure  II- 1.  (The  shaded  area  of  Figure  II- 1 shows  the  area  in 
which  events  used  in  this  evaluation  occurred.  ) During  the  contract  period 
covered  by  this  report,  data  were  available  from  twelve  of  these  stations: 
Albuquerque,  New  Mexico  (ANMO);  Mashhad,  Iran  (MAIO);  Guam,  Marianas 
Islands  (GUMO);  Narrogin,  Western  Australia  (NWAO);  Taipei,  Taiwan 
(TATO);  Wellington  (South  Karori),  New  Zealand  (SNZO);  Chiang  Mai,  Thai- 
land (CHTO);  Charters  Towers,  Australia  (CTAO);  Kabul,  Afghanistan 
(KAAO);  Matsushiro,  Japan  (MAJO);  Zongo,  Bolivia  (ZOBO);  and  Bogota, 
Columbia  (BOCO).  Stations  included  in  this  evaluation  are:  ANMO,  MAIO, 
CHTO,  CTAO,  ZOBO,  KAAO,  and  MAJO. 

Figure  II- 2 shows  diagrammatically  the  availability  of  data  for 
each  evaluated  station  from  the  start-up  date  through  March  1978.  The  ver- 
tical dashed  lines  in  Figure  II-2  enclose  the  time  frame  of  the  event  detection 
capability  studies.  Data  availability  is  nearly  continuous  for  all  stations  with- 
in their  evaluation  time  frames  with  the  exception  of  MAJO,  which  was  down 
for  extended  periods  of  time  due  to  equipment  failure. 

B.  FORMATION  OF  THE  EVENT  DATA  BASE 

This  report  presents  detection  and  discrimination  studies  for 
five  new  stations:  CHTO,  CTAO,  ZOBO,  KAAO,  and  MAJO. 
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LOCATIONS  OF  SRO  STATIONS 


AVAILABLE  SRO  DATA  FOR  EVALUATED  STATIONS  FROM 
JANUARY  1976  THROUGH  MARCH  1978 


Evaluation  of  these  stations  begins  on  station  date  1 September 
1977.  This  is  the  approximate  two  month  operational  anniversary  of  MAJO, 
which  is  the  Matsushiro,  Japan  ASRO  station,  and  the  newest  site  evaluated. 
This  two  month  waiting  period  was  initiated  in  order  to  allow  personnel  time 
to  correct  any  malfunctions  which  may  occur  early  in  the  operational  life  of 
a seismic  station.  As  of  the  middle  of  August  1978,  microfiche  records  for 
the  stations  of  interest  were  unavailable  past  station  date  31  November  1977. 
This  date  marks  the  end  of  the  new  station  evaluation  time  frame. 

Events  within  this  time  frame  were  edited  from  the  NORSAR 
and  NEIS  lists.  Those  events  satisfying  the  previously  imposed  geographic 
restraints  formed  a final  Eurasian  data  base  of  247  events. 

The  short-period  data  base  used  in  this  evaluation  is  simply 
a subset  of  the  above- described  data  base.  All  members  of  this  subset  have 
station- epicenter  separations  of  less  than  103  degrees.  This  restriction  was 
placed  on  the  short-period  data  base  since  detection  of  the  P wave  beyond  103° 
from  the  event  epicenter  is  limited  by  a shadow  zone.  At  ZOBO  this  restric- 
tion severely  limited  the  size  of  the  data  base,  since  these  stations  are  located 
far  from  the  general  area  of  interest.  Table  H-2  presents  the  mean  earthquake 
epicentral  distance  to  each  station.  (These  values  are  derived  from  the  station- 
epicenter  separation  of  the  events.  ) 

The  description  of  the  noise  sample  data  base  is  given  in  detail 
in  Section  IV  of  this  report.  In  brief,  the  short-period  noise  samples  were 
selected  from  time  gates  immediately  preceding  the  observed  signal  which 
triggered  the  automatic  detector.  Long-period  noise  samples  were  arbitrarily 
processed  at  noon  of  each  day.  Short-period  and  long-period  noise  samples 
w^re  edited  every  fourth  field  tape  day.  Noise  samples  were  quality  checked 
by  visual  examination  of  paper  plots  and/or  microfiche  and  samples  contain- 
ing signals  were  rejected. 
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c. 


DATA  PROCESSING 


Long-period  signal  and  noise  data  and  short-period  noise  data 
were  processed  in  two  stages;  a pre-analysis  processing  and  an  analysis  pro- 
cessing stage  involves  the  use  of  the  multi-purpose  program  TISSPROG 
(Schmidt,  1978).  A description  of  this  program,  as  it  processes  SRO  data, 
follows  (Figure  II- 3). 


Given  input  data  consisting  of  epicentral  locations  and  origin 
times,  TISSPROG  will  estimate  short-  or  long-period  arrival  times  and  edit 
events  or  noise  samples  from  a field  tape.  Short-period  data  are  resampled 
to  a one-tenth  of  a second  time  interval,  and  long-period  data  are  resampled 
from  a one- second  to  a two- second  time  interval.  Long-period  edit  gates  are 
automatically  set  at  4096  seconds,  and  short-period  edit  gates  are  determined 
by  the  'on  time'  of  the  short-period  detector  (see  Section  HI).  Short-period 
edit  gates  are  limited  to  204.  8 seconds.  Trace  means  are  next  removed,  and 
long-period  data  are  rotated  from  their  vertical,  north,  east  configuration  to 
a vertical,  transverse,  radial  configuration.  At  this  stage,  samples  are  sav- 
ed on  an  event  tape  for  further  analysis.  Finally,  TISSPROG  produces  0.  5- 
4.  0 Hz  bandpass  filtered  short-period  or  0.  023-0.  059  bandpass  filtered  long- 
period  plots. 

Noise  and  signal  analyses  follow.  Since  procedures  are  detail- 
ed in  later  sections,  only  brief  descriptions  will  be  included  here. 

All  data  samples  were  visually  quality  checked.  Long-period 

events  were  analyzed  for  detection  and  25  second  M - m,  discrimination  ca- 

s b 

pability.  Long-period  noise  samples  yielded  512  point  noise  analysis  gates 
These  noise  gates  were  further  processed  to  produce  peak  20,  25,  and  30  sec- 
ond noise  amplitudes,  RMS  noi3e  amplitudes  in  the  10-25,  17-41,  and  40-64 
second  spectral  bandwidths,  and  power  spectra.  Values  were  later  grouped 
as  averages  and/or  functions  of  time. 
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Since  short-period  data  are  recorded  only  when  the  SRO  detec- 
tor is  signal  activated  (Section  III),  acceptable  extended  noise  gates  are  non- 
existant.  However,  each  time  the  detector  triggers,  the  preceding  20  seconds 
of  data  are  recorded  as  a lead-in  buffer.  Therefore,  the  first  12.  8 seconds 
of  short-period  detections  were  used  for  noise  analysis.  Following  visual  in- 
spection, an  analysis  processing  routine  was  employed  to  produce  a peak  one 
second  amplitude,  0.  5-4.  0 Hz  RMS  amplitude,  and  power  spectrum.  These 
values  were  also  grouped  as  averages  and/or  time  functions. 

Short- period  detection  capability  was  estimated  with  the  aid  of 
microfiche  analyses.  Computer  processing  proved  necessary  only  for  those 
few  events  satisfying  the  following  criteria: 

• No  signal  could  be  identified  visually  on  microfiche  records. 

• Estimated  signal  arrival  time  occurred  more  than  30  seconds 
after  the  start  time  of  the  first  record  of  an  automatic  detec- 
tion. 

D.  PROCESSING  SUMMARY 

Table  II-3  summarizes  the  processing  results  of  the  work  per- 
formed during  the  current  contract  period  on  the  Seismic  Research  Observa- 
tories evaluation  task.  In  this  table,  the  'SP'  and  'LP'  under  the  heading 
'DATA  TYPE'  refer  to  short-period  data  and  long-period  data,  respectively. 
The  heading  'DETECTED'  refers  to  the  number  of  events  which  were  visually 
detected  (Calcomp  plot  for  LP,  microfiche  or  Calcomp  plot  for  SP)  under  the 
detection  criteria  of  Section  V.  The  heading  'NOT  DETECTED'  refers  to  the 
number  of  events  for  which  only  seismic  noise  was  observed  in  the  signal  gate. 
The  heading  'MIXED'  refers  to  the  number  of  events  which  were  not  detectable 
due  to  the  presence  of  some  other  signal  in  the  signal  gate.  The  heading 
'MICROFICHE  DETECTED  ONLY'  refers  to  the  number  of  events  not  detected 
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by  the  short-period  automatic  detector.  The  heading  'SYSTEM  FAILURE  - 
NO  DATA  RECORDED'  refers  to  the  number  of  events  for  which  no  data  were 
recorded  since  the  system  had  gone  down  (failed).  The  heading  'MALFUNC- 
TION' refers  to  the  number  of  events  for  which  the  detection  status  could  not 
be  determined  due  to  malfunctions  (spikes,  glitches,  and  the  like)  in  the  sig- 
nal gate.  The  heading  'A  < 103°'  refers  to  the  number  of  events  for  which 
short-period  data  were  not  processed  due  to  their  epicentral  distances  being 
greater  than  103°. 

The  quality  of  the  data  recorded  at  the  Seismic  Research  Ob- 
servatories was  very  good  in  most  cases.  The  types  of  malfunctions  which 
degrade  the  data  quality  were  described  in  detail  in  the  preliminary  evalua- 
tion of  the  Seismic  Research  Observatories  (Strauss,  1976).  The  overall  fre- 
quency of  occurrence  of  these  malfunctions  can  be  estimated  from  Table  II- 3. 

Long-period  malfunctions  which  occurred  at  ZOBO  and  MAJO 
were  hardware  related,  and  these  stations  have  since  been  repaired. 

Estimates  of  station  reliability  were  made  using  the  following 
argument.  If  a station  is  perfect  in  that  it  always  produces  seismic  record- 
ings which  the  analyst  can  check  for  detections,  it  is  considered  to  have  a 
reliability  factor  of  1. 0.  If  on  the  other  hand,  the  station  is  never  producing 
seismic  recordings  due  to  instrument  problems  (instruments  recording  im- 
properly or  not  at  all),  it  is  considered  to  have  a reliability  factor  of  0.  0. 

In  practice,  of  course,  the  reliability  factor  lies  somewhere  between  these 
extremes.  Since  station  down  time  and  station  malfunction  time  are  the  two 
factors  which  render  the  station  reliability  less  than  1. 0,  tne  reliability  fac- 
tor is  defined  as 

Reliability  Factor  = 1.0-  (Percentage  of  time  station  is  down 

+ Percentage  of  time  station  is  mal- 
functioning) 
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where  the  two  percentages  are  estimated  empirically  from  the  data  of  Table 
II- 3.  Thus,  the  percentage  of  time  the  station  is  down  is  estimated  from  the 
number  of  events  for  which  no  data  were  recorded  divided  by  the  total  number 
of  events  for  which  processing  was  attempted  and  the  percentage  of  time  the 
station  is  malfunctioning  is  the  number  of  events  for  which  malfunctions 
(spikes,  glitches,  and  data  drop-outs)  overrode  the  seismic  data  divided  by 
the  total  number  of  events  for  which  processing  was  attempted.  When  com- 
puting the  short-period  reliability  factors,  those  events  for  which  no  data 
were  recorded  due  to  the  automatic  detector  were  excluded,  since  it  is  not 
possible  to  determine  how  many  would  contain  malfunctions. 

The  reliability  estimates  are  presented  in  Table  II-4.  The 
relatively  low  reliabilities  of  MAJO  and  ZOBO  are  due  primarily  to  station 
down  time,  when  no  data  were  recorded  (see  Figure  H-2). 

From  the  data  base  processed  in  this  evaluation,  it  was  pos- 
sible to  form  estimates  of  the  probability  than  an  event  of  interest  will  be 
masked  by  some  other  event  (a  'mixed'  event).  From  the  data  of  Table  II- 3, 
one  can  see  that  the  probability  that  any  short-period  event  will  be  mixed  is 
very  small.  This  low  probability  of  mixing  for  short-period  data  is  due  to 
the  length  of  the  signal  coda,  which  is  quite  short  relative  to  the  average  time 
between  signals. 

The  long-period  data  presents  quite  a different  picture.  Due 
to  their  relatively  long  signal  coda,  there  is  a significant  probability  of  mix- 
ing for  long-period  data.  Table  II-5  presents  estimates  of  the  probability  of 
mixing  derived  from  the  processed  data  base.  The  high  frequency  of  mixed 
events  evident  at  ZOBO  is  most  probably  due  to  the  distance  from  that  station 
to  Eurasia.  The  signal  arrives  over  a large  gate  consequently,  there  is  more 
opportunity  for  mixing  to  occur. 
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TABLE  II- 4 

STATION  RELIABILITY  ESTIMATES 


Station 

Short- Period 

Long  - Period 

CHTO 

0.  85 

0.84 

CTAO 

0.  99 

0.  99 

ZOBO 

0.  72 

0.63 

KAAO 

0.  93 

0.89 

MAJO 

0.  71 

0.66 

SECTION  in 

THE  SHORT- PERIOD  AUTOMATIC  DETECTOR 


For  convenience  in  the  SRO  short- period  automatic  detector 
performance  evaluation,  the  detector  performance  criteria  (Subsection  A)  and 
the  detector  design  description  (Subsection  B)  from  last  year's  SRO  evaluation 
report  (Strauss  and  Weltman,  1977)  are  repeated.  These  subsections  are  fol- 
lowed by  the  current  evaluation  of  the  detector  performance. 


A.  DISCUSSION 


This  section  reviews  the  design  of  the  short-period  automatic 
detector  currently  in  operation  at  the  Seismic  Research  Observatory  stations 
and  discusses  the  problems  of  this  detector  discovered  during  the  course  of  this 
evaluation. 


The  prime  reason  for  incorporating  a short-period  automatic 
detector  at  the  Seismic  Research  Observatory  stations  was  to  minimize  the 
amount  of  short-period  data  recorded  during  which  no  short-period  seismic 
events  could  be  observed  and  thus  maximize  the  recording  life  of  each  mag- 
netic tape.  This  would  minimize  expenditures  for  the  purchasing  and  trans- 
portation of  magnetic  tapes  and  the  amount  of  data  lost  during  tape  changes 
without  seriously  decreasing  the  number  of  short-period  seismic  events  re- 
corded. Since  the  single  component  of  short-period  motion  is  sampled  twenty 
times  per  second  while  the  three  components  of  long-period  motion  are  each 
sampled  once  per  second,  the  bulk  of  each  tape  (87%)  would  be  filled  with 
short-period  data  if  short-period  data  were  recorded  continuously.  In  this 
continuous  recording  mode,  the  recording  life  of  a 2400  foot,  800  BPI  mag- 
netic tape  would  be  approximately  four  days.  With  the  automatic  detector  in 
operation,  this  has  been  extended  to  up  to  sixteen  days. 
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The  design  of  any  automatic  detector  should  take  into  account 
the  following  points: 

• Maximum  detection  probability  consistent  with  the  largest 
tolerable  false  alarm  rate 

• Stable  false  alarm  rate  over  a wide  range  of  ambient  noise 
levels 

• Minimum  number  of  detections  declared  from  the  same  event 
while  retaining  the  ability  to  detect  possible  secondary  phase 
arrivals  within  the  coda 

• Minimum  time  necessary  to  re-establish  equilibrium  after 
parameter  value  changes 

• Low  sensitivity  to  degraded  data  quality  (Swindell  and  Snell, 
1977). 

In  this  evaluation  of  the  automatic  detector  design,  data  output 
by  the  automatic  detector  will  be  assessed  in  terms  of  the  above  points. 

B.  DESIGN  OF  THE  PRESENT  SHORT- PERIOD  AUTOMATIC  DETECTOR 

Each  Seismic  Research  Observatory  station  is  equipped  with  a 
short-period  automatic  event  detector  (hereafter  referred  to  as  the  SRO  auto- 
matic detector)  which  is  intended  to  discriminate  between  seismic  noise  and 
short-period  seismic  events.  The  design  of  the  SRO  automatic  detector  is 
bas  ’d  on  a study  carried  out  by  the  Charles  Stark  Draper  Laboratory  (Eterno, 
et  al.  , 1974). 

The  design  concept  of  the  SRO  automatic  detector  is  illustrated 
by  Figure  III- 1 . The  short- period  analog  data  from  the  sensor  unit  is  first 
digitized  at  a rate  of  twenty  samples  per  second.  At  the  SRO  automatic  detector 
input,  each  datum  is  converted  to  a twelve  bit  integer  quantity.  The  digitized 
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FIGURE  EI-l 

SHORT- PERIOD  DETECTOR  CONCEPT 
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data  next  passes  through  a band-limiting  four  pole-pair  Butterworth  filter 
with  two  zeros  at  zero  frequency.  (This  filter  is  constructed  by  cascading 
four  re-.utJ.vfc  single  pole- pair  filters.  A single-zero  filter  is  placed  in  front 
of  the  first  and  fourth  cascaded  filters.  ) The  lower  cutoff  of  this  filter  is  0.  5 
Hz,  the  upper  cutoff  is  2.  0 Hz,  the  filter  response  within  the  passband  is  flat 
to  within  ± 1 dB,  and  the  rejection  rate  outside  the  passband  is  more  than  25 
dB  /octave. 

The  data  are  then  rectified  and  passed  to  the  short-term  aver- 
ager (STA).  Utilizing  a single-pole  recursive  stage,  one- second  exponentially- 
weighted  short-term  averages  are  constantly  generated.  The  short-term  aver- 
ager imparts  an  absolute  gain  of  approximately  13  dB  to  the  data.  The  short- 
term averages  are  combined  to  form  a long-term  average  (LTA).  Long-term 
averages  may  be  formed  over  operator- selected  time  intervals  of  10,  20,  30, 
or  40  seconds.  The  approximate  13  dB  gain  imparted  by  the  short-term  aver- 
ager allows  greater  accuracy  in  calculating  the  long-term  average.  However, 
once  the  calculation  is  complete,  the  long-term  average  is  re-scaled  to  unitize 
the  gain  at  zero  frequency. 

The  SRO  automatic  detector  tests  for  detections  by  continuous- 
ly comparing  the  short-term  average  with  the  long-term  average.  The  long- 
term average  is  first  scaled  upward  by  an  operator- selected  threshold  con- 
stant (the  turn-on'  threshold).  This  scaling  is  intended  to  minimize  false 
alarms.  If  the  short-term  average  rises  above  the  scaled  long-term  average, 
the  long-term  average  is  'frozen'  at  its  current' value.  If  the  short-term  aver- 
age than  remains  above  the  frozen  value  of  the  long-term  average  for  one  sec- 
ond (twenty  samples)  a detection  is  declared  and  the  short-period  data,  which 
hasv  been  held  in  a separate  memory  buffer,  is  written  on  the  magnetic  tape, 
a,  continuous  recording  begins.  If  the  short-term  average  drops  below  the 
frozen  long-term  average  before  the  one  second  gate  is  over,  no  detection  is 
de<  lared,  r.he  long-term  average  is  'unfrozen'  and  updated. 


IH-4 


Once  a detection  has  been  declared,  the  frozen  long-term  average  is  scaled 
upward  by  a second  operator- selected  threshold  constant  (the  'turn-off'  thresh- 
old). When  the  short-term  average  drops  below  the  rescaled  long-term  aver- 
age, the  event  is  considered  to  be  over  and  the  current  buffer  of  short-period 
data  completes  the  magnetic  tape  record.  Scaling  the  frozen  long-term  aver- 
age by  the  turn-off  threshold  is  intended  to  minimize  the  recording  of  the 
sometimes  persistent  coda  that  can  follow  the  initial  arrival  of  a short-period 
event. 

Thus,  the  operator  can  ’tune'  the  SRO  automatic  detector  via 
four  adjustable  constants:  the  long-term  averager  time  constant,  the  detection 
turn-on  threshold,  the  detection  turn-off  threshold,  and  the  digital  gain  factor. 
In  addition,  the  operator  can  override  the  SRO  automatic  detector  and  record 
short-period  continuously. 

Set  procedures  exist  for  setting  the  proper  digital  gain  factor. 
However,  selection  of  the  other  constants  is  largely  a matter  of  trial-and- 
error.  The  operator  zeros  in  on  optimum  values  after  selecting  initial  values 
within  experimentally  determined  ranges. 

C.  EVALUATION  OF  THE  AUTOMATIC  DETECTOR 

Microfiches,  covering  the  time  period  1 September  1977 
through  30  November  1977,  were  obtained  for  SRO  station  CHTO  (Chiang  Mai, 
Thailand).  Every  third  day  was  visually  examined  and  analyst  detection  logs 
were  compiled.  Detection  criteria  were  a 3 dB  peak- signal-to-peak  noise 
ratio  and  an  apparent  signal  character  as  determined  by  the  analyst.  These 
logs  were  compared  to  automatic  detector  logs  which  were  generated  from 
CHTO  field  tapes. 

The  automatic  detector  logs  are  based  on  the  conformity  of 
record  times  rather  than  on  the  detection  status  bit  included  in  each  field 
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tape  short-period  record.  This  method  was  preferred  since  detections  of 
less  than  one  record  duration  (49  seconds)  may  not  be  flagged  by  the  status 
bit. 

Comparison  of  the  analyst- gene  rated  and  automatic  detector 
log  yielded  detection  statistics  which  may  be  used  to  gauge  detector  perfor- 
mance. Semi-monthly  sums  over  every  third  day  (a  total  of  5 days  or  120 
hours  per  period)  of  automatic  detector  'misses',  'false  alarms',  and  'con- 
firmed detections'  are  given  in  Table  III- 1 . An  item  in  the  column  labeled 
'Number  Detected  by  Analyst  Only'  is  called  a 'miss'  and  refers  to  the  condi- 
tion: signal  in  the  appropriate  time  gate,  no  detection  declared.  The  column 
labeled  'Number  Detected  by  Automatic  Detector'  reflects  the  false  alarm 
rate  and  refers  to  the  condition:  no  signal  in  the  time  gate,  detection  declared. 
The  column  labeled  'Number  Detected  by  Both'  shows  the  number  of  automatic 
detections  confirmed  by  the  analyst.  The  total  number  of  valid  detections  then 
is  found  by  adding  the  numbers  of  missed  and  confirmed  detections.  The  num- 
bers in  Table  ni-1  translate  into  the  detector  operating  characteristics.  Fig- 
ure III- 2. 

We  observe  that  the  probability  of  detection,  P^,  ranges  ap- 
proximately from  50%  to  80%  at  false  alarm  rates  (FAR)  varying  between  0.  2 
and  2.  1 false  alarms  per  hour  (FA/H).  While  the  averaging  periods  (5  days 
in  a half- month  interval)  should  be  sufficiently  long  to  give  solid  statistics, 
the  probability  of  detection  seems  rather  independent  of  the  FAR.  This  prob- 
ability reflects  the  FAR  instability  inherent  in  a straight  STA/L.TA  detector 
which  does  not  account  for  changes  in  the  STA  variance  for  noise  (Lacoss, 

1972;  Swindell  and  Snell,  1977;  Unger,  1978). 

Several  points  should  be  considered  with  regard  to  the  validity 
of  the  analyst- detector  comparisons.  The  analyst  was  not  concerned  with 
source -to- station  distances  when  picking  signals.  Since  the  automatic  detec- 
tor is  not  tuned  for  local  events  (Personal  Communication,  John  Hoffman, 
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1978),  the  number  of  missed  events  may  be  higher  than  what  would  be  expect- 
ed for  teleseismic  data  only.  Also,  as  illustrated  in  Figure  III- 1 , a bandpass 
filter  precedes  the  automatic  detection  process.  The  microfiche  viewed  was 
not  bandpass  filtered,  and  the  large  number  of  analyst  picks,  along  with  limit- 
ed computer  availability,  prohibited  such  processing.  Thus,  the  estimated 
operating  characteristics  may  be  biased.  Finally,  as  can  be  seen  in  Table 
III-l,  the  CHTO  automatic  detector  malfunctioned  for  eight  days.  Data  were 
recorded  almost  continuously  and  no  attempt  was  made  at  automatic  detector 
evaluation  for  that  time  period. 

Another  problem  causing  FAR  instability,  encountered  in  last 
year's  evaluation  (Strauss  and  Weltman,  1977)  is  illustrated  in  Table  IH-2. 

In  many  instances,  there  was  no  relation  between  the  number  of  detections 
and  the  amount  of  recording  time,  which  indicated  possible  difficulties  in  the 
setting  and  freezing  of  detector  parameters  (Strauss  and  Weltman,  1977). 

An  identical  study  was  conducted  this  year  on  the  five  new  sta- 
tions  under  evaluation.  The  percentage  time  recorded  and  the  number  of  de- 
tections were  charted  for  a nine-day  period  and  are  shown  in  Table  III- 3.  The 
data  seem  much  more  constant  with  time  for  a given  station.  This  suggests 
that  an  improvement  has  been  made  in  the  parameter  setting  and  freezing 
problem. 
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TABLE  HI- 2 


PERFORMANCE  OF  AUTOMATIC  DETECTOR  FOR  AN  ARBITRARY 

TEN- DAY  PERIOD  IN  1976 
(Strauss  and  Weltman,  1977) 


Date 

Number  of  Detections 

ANMO 

GUMO 

MAIO 

NWAO 

TATO 

SNZO 

06/21/76 

6 

D 

11 

23 

9 

7 

06/22/76 

19 

9 

0 

1 

1 

06/23/76 

14 

n 

9 

3 

4 

1 

06/24/76 

16 

6 

1 

2 

2 

06/25/76 

16 

Bf 

20 

6 

0 

2 

06/26/76 

4 

8 

15 

3 

0 

0 

06/27/76 

8 

18 

11 

4 

2 

3 

06/28/76 

15 

64 

12 

3 

1 

2 

06/29/76 

13 

„ 62 

6 

4 

2 

6 

06/30/76 

15 

22 

5 

3 

7 

4 

Date 

Percentage  of  Time  Recorded 

ANMO 

GUMO 

MAIO 

NWAO 

TATO 

SNZO 

06/21/76 

59 

10 

12 

1 

2 

D 

06/22/76 

6 

33 

0 

9 

E9 

06/23/76 

3 

42 

>1 

1 

>i 

06/24/76 

4 

0 

>1 

>1 

>i 

06/25/76 

6 

3 

14 

20 

0 

89 

06/26/76 

1 

1 

32 

3 

0 

100 

06/27/76 

2 

2 

26 

2 

>1 

64 

06/28/76 

3 

9 

22 

1 

96 

>1 

06/29/76 

3 

49 

11 

>1 

5 

42 

06/30/76 

3 

2 

54 

>1 

1 

41 
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TABLE  in- 3 


PERFORMANCE  OF  AUTOMATIC  DETECTOR  FOR 
AN  ARBITRARY  NINE-DAY  PERIOD  IN  1977 


Date 

Number  of  Detections 

CHTO 

CTAO 

KAAO 

MAJO 

ZOBO 

11/14/77 

13 

20 

28 

8 

26 

11/15/77 

4 

29 

30 

6 

22 

11/16/77 

19 

24 

41 

10 

37 

11/17/77 

12 

22 

24 

2 

33 

11/18/77 

15 

20 

33 

11 

29 

11/19/77 

13 

18 

31 

4 

33 

11/20/77 

7 

56 

29 

7 

35 

11/21/77 

10 

42 

32 

14 

50 

11/22/77 

10 

28 

38 

4 

26 

Date 

Percentage  of  Time  Recorded 

CHTO 

CTAO 

KAAO 

MAJO 

ZOBO 

11/14/77 

1 

3 

10 

2 

6 

11/15/77 

1 

3 

7 

3 

2 

11/16/77 

11 

4 

7 

5 

6 

11/17/77 

3 

3 

4 

2 

5 

11/18/77 

8 

3 

11 

2 

5 

11/19/77 

1 

2 

7 

>1 

7 

11/20/77 

1 

5 

6 

1 

4 

11/21/77 

4 

4 

11 

4 

9 

11/22/77 

4 

7 

11 

1 

4 
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SECTION  IV 
NOISE  ANALYSIS 


A.  DISCUSSION 

The  goal  of  this  section  is  to  characterize  the  noise  field  at 
each  site  under  evaluation.  Presented  in  this  section  are  peak  1,  20,  25, 
and  30  second  noise  values,  RMS  trends  and  average  RMS  values  in  the  0.  5- 
1.25,  10-25,  17-41,  and  40-64  second  passbands,  and  average  short-  and 
long-period  RMS  spectra. 

All  noise  values  in  this  section  are  presented  without  instru- 
ment response  correction,  since  it  is  recognized  that  the  analyst  is  primarily 
concerned  with  the  noise  as  he  will  see  it,  i.  e.  , after  it  has  passed  through 
the  sensing,  filtering,  and  recording  instrumentation.  Instrument  response 
corrected  noise  values  are  presented  in  Appendix  B. 

B.  VERTICAL  COMPONENT  SHORT- PERIOD  NOISE 

The  short-period  noise  analysis  for  the  five  new  stations  cur- 
rently being  evaluated  (CHTO,  CTAO,  KAAO,  MAJO,  and  ZOBO)  extended 
from  1 July  1977  to  31  December  1977.  Subject  to  data  availability,  three 
samples  were  taken  from  every  fourth  day,  and-from  these  three  samples, 
one  was  chosen  to  represent  that  day.  All  samples  were  picked  in  one  of  two 
ways.  First,  on  days  that  contained  events  from  the  data  base,  samples 
were  taken  from  the  portion  of  the  automatic  detector  edit  preceeding  a de- 
tected event.  On. days  that  contained  no  events  from  the  data  base,  samples 
were  picked  randomly  from  automatic  detector  logs  (see  Section  III)  with  the 
only  requirement  being  that  the  time  separating  consecutive  automatic  detec- 
tions exceed  one  hour.  This  was  done  to  prevent  selecting  a sample  inside 
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a series  of  multiple  detections.  The  standardized  length  of  each  noise  sam- 
ple was  12.  8 seconds,  starting  on  the  first  point  of  the  detector  edit. 


After  all  unacceptable  samples  were  screened  out,  each  sam- 
ple was  filtered  by  a 0.  5-4.  0 Hz  bandpass  filter.  RMS  noise  was  then  com- 
puted by  the  equation  i 


RMS  NOISE 


r n 2 t 

£ <*/ 
i=l 

. n 


where 


n = number  of  data  points 


x.  = the  i data  point. 

l 


The  calculated  RMS  values  (in  millimicrons)  were  then  plotted  against  Julian 
day  in  Figures  IV- 1 to  IV- 3 and  the  monthly  RMS  noise  trends  derived  from 
these  values  were  plotted  in  Figures  IV-4  to  IV-6. 

The  RMS  noise  trends  for  stations  CHTO  and  ZOBO  showed  little 
change  with  time.  Of  these  two,  ZOBO  had  the  lowest  and  the  most  consistent 
noise  level.  CHTO,  on  the  other  hand,  showed  several  small  peaks  between 
days  200  and  270.  These  are  possibly  due  to  periods  of  storm  activity. 

Station  KAAO  showed  very  good  day-to-day  consistency,  but 
the  trend  plot  showed  a peak  through  August  and  September.  Stations  CTAO 
and  MAJO  showed  the  most  variability  in  both  daily  and  trend  RMS  noise, 
wit*’  MAJO  increasing  slowly  with  time  and  CTAO  decreasing  with  time. 

Table  IV-1  contains  the  mean  short-period  RMS  noise  values. 

A correlation  between  mean  RMS  noise  level  at  a station  and  the  distance  be- 
tw-  ,^n  that  station  and  the  nearest  coast  is  evident.  This  correlation  was  pre- 
viously noted  by  Strauss  and  Weltman  (1977).  The  stations  with  the  lowest 
noise  levels  are  located  farthest  from  the  nearest  coast, 


l 


IV- 2 


CHTO  AND  CTAO  SHORT- PERIOD  RMS  NOISE 


ZOBO  SITE  51 


-PERIOD  RMS  NOISE 


SITE  33 
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FIGURE  IV- 4 

CHTO  AND  CTAO  SHORT- PERIOD  RMS  NOISE  TRENDS 
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FIGURE  IV-6 

MAJO  SHORT- PERIOD  RMS  NOISE  TRENDS 
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MEAN  SHORT- PERIOD  RMS  NOISE  (VERTICAL  COMPONENT) 


Station 

Approximate  Distance 
to  Nearest  Coast  (km) 

Mean  RMS  Noise 
(m/i) 

Standard 

Deviation 

Number  of 
Sample  s 

CHTO 

290 

1. 67 

0.  67 

34 

CTAO 

120 

5.  55 

2.  2Z 

28 

KAAO 

1100 

1.  94 

0.  69 

37 

MAJO 

110 

3.  53 

1.  31 

28 

ZOBO 

300 

1. 02 

0.  36 

37 

i.  e.  , ZOBO,  CHTO,  and  KAAO.  MAJO  and  CTAO  are  much  closer  to  the 
coast  than  the  other  three  stations  and  their  noise  levels  are  significantly 
higher.  As  can  be  seen  from  Table  IV- 1,  however,  the  correlation  is  not 
exact.  For  this  reason,  the  'baseline  RMS',  the  level  above  which  short- 
period  noise  is  caused  by  sea  wave  energy,  cannot  be  accurately  defined. 
These  data  place  the  baseline  between  1.  0 m/x  and  2.  0 m/x  . 

Maximum  zero-to-peak  one- second  period  noise  was  also  mea- 
sured for  each  sample.  The  statistics  of  these  measurements  are  presented 
in  Table  IV-2  in  terms  of  the  mean  and  standard  deviation  of  the  logarithm  of 
the  measured  values. 

Short-period  noise  analysis  was  expanded  this  year  to  include 
spectral  analysis.  Each  sample  was  filtered  with  a 0.  5-4.  0 Hz  bandpass  fil- 
ter and  the  mean  and  standard  deviation  of  the  amplitude  spectrum  was  cal- 
culated for  each  frequency  increment.  The  mean  spectral  density  and  the 
logarithm  of  the  mean  spectral  density  with  standard  deviations  are  plotted 
in  Figures  IV- 7 to  IV - 1 1 . 

In  a study  of  ambient  earth  motion.  Fix  (1972)  observed  peaks 
in  the  spectral  density  at  periods  of  0.  3 and  0.  5 seconds,  with  the  0.  5 second 
peak  being  slightly  less  prominent.  These  were  proved  to  be  primarily 
fundamental  mode  Rayleigh  waves  (Douze,  1967).  The  results  of  our  analyses 
compare  fairly  well  with  those  findings.  MAJO  and  CHTO  show  peaks  at  ap- 
proximately 0.  3 seconds,  and  ZOBO  and  CHTO  show  clear  peaks  at  0.  5 and 
0.  55  seconds,  respectively.  Stations  CTAO  and  KAAO,  however,  show 
neither  of  these  peaks  clearly.  KAAO  shows  a peak  at  approximately  0.  75 
seconds,  but  at  shorter  periods  the  spectrum  contains  no  clear  peaks.  There 
is  a concentration  of  energy  between  periods  of  0.  3 to  0.  6 seconds,  but  there 
is  ,iO  'structure'  to  the  spectrum  like  those  at  the  other  stations  mentioned. 
This  concentration  is  probably  caused  by  cultural  noise.  The  CTAO  spectrum 
is  ever,  more  unlike  the  spectra  of  the  other  stations.  Not  only  are  the 
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TABLE  IV- Z 


SHORT- PERIOD  NOISE 

LOGj Q (PEAK  ONE-SECOND  NOISE  AMPLITUDE)  STATISTICS 


Station 

Mean  Logj  q 
Peak  1- Second 

Standard  Deviation 
LogjQ  Peak  1- Second 

Number  of 
Samples 

CHTO 

0.  369 

0.  182 

29 

CTAO 

0.  614 

0.  260 

21 

• 

KAAO 

0.463 

0.  207 

33 

MAJO 

0.  556 

0.  213 

26 

ZOBO 

0.  098 

0.  199 

33 
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SITE  33 


AVERAGE  RMS  AMPLITUDE  SPECTRA  CHTO  SHORT-PERIOD  NOISE 
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AVERAGE  RMS  AMPLITUDE  SPECTRA  CTAO  SHORT- PERIOD  NOISE 
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SITE  52 
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AVERAGE  RMS  AMPLITUDE  SPECTRA  KAAO  SHORT -PERIOD  NOISE 


SITE  53 
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amplitudes  inflated,  but  there  are  no  apparent  concentrations  of  energy  at 
any  frequency.  No  ready  explanation  can  be  given  for  this  peculiarity.  Per- 
haps the  site  is  contaminated  by  some  strong  source  of  cultural  noise, 

C.  THREE  COMPONENT  LONG- PERIOD  NOISE 

The  goals  of  the  long-period  noise  analysis  were  to  estimate 
long-period  RMS  noise  levels,  peak  noise  amplitudes,  and  the  spectral  con- 
tent of  the  noise  field  for  each  of  the  three  components  (V,  N,  E)  at  each  sta- 
tion under  evaluation. 

Long-period  noise  processing  usually  began  two  months  after 

each  station's  declared  operational  date  and  continued  as  long  as  time  and 

computer  availability  permitted.  The  time  periods  covered  are  shown  in 

Table  IV- 3.  The  two-months  delay  was  instituted  to  ensure  that  each  station 

was  functioning  properly  prior  to  evaluation.  Noise  analysis  time  frames 

were  limited  to  one  data  year  and  periods  of  known  hardware  malfunction 

% 

were  deleted  from  the  noise  data  base. 

The  4096- second  noise  samples  were  processed  as  described 
in  Section  II.  The  V,  N,  E configuration  was  maintained.  Samples  were 
processed  at  1200  hours  every  fourth  day,  where  possible.  Samples  were 
visually  screened  for  signals  and  unreported  system  malfunctions  and  1024- 
second  noise  analysis  gates  were  selected.  If  an  entire  sample  proved  un- 
acceptable, a second  attempt  followed,  and  an  acceptable  noise  sample  was 
usually  found  at  an  edit  time  within  two  hours  of  the  first.  Noise  samples 
and  analysis  gates  were  then  input  into  a primary  analysis  program  which 
performed  the  following  functions: 

• Computed  RMS  noise  values  with  and  without  instrument  response 

in  the  following  passbands:  10-25,  17-41,  and  40-64  seconds. 

• Measured  zero-to-peak  20,  25,  and  30- second  noise  amplitudes. 
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TABLE  IV- 3 

LONG- PERIOD  NOISE  DATA  BASE 


Station 

Start  Date 
(Year- Julian  Day) 

End  Date 

(Year- Julian  Day) 

Number  of 
Samples 

AN  MO 

77009 

76362 

88 

CHTO 

77251 

78022 

s 31 

MAIO 

76001 

76365 

82 

CTAO 

77003 

77363 

83 

ZOBO 

77153 

77341 

39 

KAaG 

77191 

78050 

53 

MAJO 

78090 

78177 

20 
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Computed  power  spectra. 


The  variety  of  RMS  and  peak  amplitude  measurements  were  made  for  the 
following  reasons.  The  17-41  second  passband  encompasses  the  signal  band 
traditionally  of  interest  from  a detection  and  discrimination  capability  view- 
point; the  10-25  second  passband  is  centered  upon  the  17-second  microseis- 
mic  peak  evident  at  most  stations,  and  the  40-64  second  passband  provides 
a look  at  lower  frequency  noise  present  at  SRO  sites.  The  variety  of  peak 
amplitude  measurements  further  describes  the  noise  field  at  each  site. 

All  RMS  and  peak  noise  values,  and  the  first  half  of  each  power 
spectrum  (128  points)  were  recorded  on  magnetic  tape.  This  tape  served  as 
input  to  the  final  analysis  and  plotting  routines. 

Long-period  RMS  noise  values  (uncorrected  for  instrument 
response)  for  the  vertical,  north,  and  east  components  are  plotted  versus 
Julian  day  in  Figures  IV-12a  through  IV-17c  for  stations  ANMO,  MAIO, 

CHTO,  CTAO,  ZOBO,  and  KAAO.  Data  gaps  of  up  to  eight  days  usually  re- 
flect an  inability  to  find  uncontaminated  noise  samples.  Larger  gaps  usually 
imply  station  down-time. 

Figures  IV-12a,  IV-12b,  and  IV-12c  show  ANMO  RMS  noise 
versus  Julian  day  within  the  three  passbands.  For  both  ANMO  and  MAIO, 
which  are  represented  by  Figures  IV-14a,  IV-14b,  and  IV-14c,  the  greater 
variances  in  daily  noise  values  appear  in  days  1 to  90  and  again  in  days  300 
to  360,  implying  winter  storm  activity  at  these  stations  (Strauss  and  Weltman, 

1977).  The  most  stable  noise  passband  for  ANMO  and  MAIO  is  the  40-64 
second  passband,  a fact  common  to  most  stations.  MAIO  exhibits  another 
common  property;  vertical  component  noise  is  significantly  lower  in  ampli- 
tude than  horizontal  component  noise  in  the  40-64  second  passband.  CHTO, 

Figures  IV- 13a,  IV- 13b,  and  IV- 13c,  shows  no  real  pattern;  its  daily  vari- 
ance remains  rather  constant.  RMS  values  seem  most  stable  in  the  40-64 
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SITE  30 
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V-20 


ANMO  10-25  SECOND  RMS  NOISE 


SITE  30 

Vertical 
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ANMO  40-64  SECOND  RMS  NOISE 


CHTO  10-25  SECOND  RMS  NOISE 


CHTO  17-41  SECOND  RMS  NOISE 
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MAIO  10-25  SECOND  RMS  NOISE 


MAIO  17-41  SECOND  RMS  NOISE 
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RMS  NOISE 


CTAO  10-25  SECOND  RMS  NOISE 


CTAO  17-41  SECOND  RMS  NOISE 
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ZOBO  10-25  SECOND  RMS  NOISE 


SITE  51 

Vertical 


ZOBO  40-64  SECOND  RMS  NOISE 


KAAO  40-64  SECOND  RMS  NOISE 


second  passband,  as  they  did  for  ANMO  and  MAIO.  CTAO  (Figures  IV-15a, 
IV-15b,  and  IV-15c)  is  almost  stable  in  the  40-64  second  passband  and  shares 
MAIO's  characteristic  low  vertical  component  noise  amplitude  in  that  pass- 
band.  The  first  60  days  of  the  year  show  relatively  little  daily  variance  com- 
pared to  the  rest,  a phenomenon  probably  related  to  seasonal  weather  trends. 
ZOBO,  represented  by  Figures  IV-16a,  IV-16b,  and  IV-16c,  shows  no  notice- 
able changes  in  noise  variance,  but  40-64  second  noise  amplitudes  are  again 
noticeably  low  on  the  vertical  component.  KAAO,  Figures  IV-17a,  IV-17b, 
and  IV-17c,  is  another  station  which  displays  little  change  in  noise  variance 
with  time.  The  low  vertical  noise  phenomenon  is  apparent  in  all  passbands 
for  this  station  and  the  40-64  second  passband  again  displays  the  least  vari- 
ance. 

Station  MAJO  suffered  from  hardware  malfunctions  for  signifi- 
cant portions  of  the  data  period  considered,  and  long-period  noise  trend  anal- 
ysis for  this  station  has  been  postponed.  Twenty  noise  samples  were  proces- 
sed from  the  later  time  period  indicated  in  Table  IV- 3.  These  samples  were 
used  to  generate  preliminary  average  values.  Mean  RMS  noise  values  and 
associated  standard  deviations  for  the  seven  stations  are  listed  in  Tables 
IV-4,  IV-5,  and  IV- 6.  These  tables  represent  the  10-25,  17-41,  and  40-64 
second  passbands,  respectively.  The  lack  of  strong  correlation  between 
long-period  noise  amplitude  and  station-to-nearest-coastline  distance  found 
last  year  is  again  apparent  (Strauss  and  Weltman,  1977).  Once  more,  the 
stations  readily  divide  into  a higher  noise  group:  ANMO,  CHTO,  and  CTAO; 
and  a lower  noise  group:  MAIO,  ZOBO,  and  KAAO.  Interestingly,  last 
year  ANMO  was  a member  of  the  low  noise  group.  At  any  rate,  CHTO  and 
CTAO  are  nearer  the  ocean  than  MAIO,  ZOBO,  and  KAAO,  so  wave- injected 
energy  is  again  implicated  as  a major  noise  factor  for  some  coastal  stations, 
allhough  it  is  certainly  not  the  only  one. 
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TABLE  IV- 4 


MEAN  10-25  SECOND  RMS  NOISE  AMPLITUDES  IN  m/i 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

11. 56 

4.  33 

10.  96 

4.  41 

11. 67 

" " 

4.  00 

CHTO 

11.16 

3.  54 

11.  70 

4.47 

12.  27 

4.  42 

MAIO 

8.  14 

3.  46 

7.  11 

2.  78 

7.  19 

% 

3.  30 

CTAO 

12.  22 

3.  89 

12.  53 

3.  47 

11.  00 

3.  74 

ZOBO 

8.  23 

2.  79 

8.  23 

2.  86 

8.  15 

2.  46 

KAAO 

8.  96 

3.  62 

9.55 

4.  06 

8.  65 

3.  38 

MAJO 

6.  64 

2.  04 

7.  68 

2.  01 

7.  79 

2.  40 

S.  D.  = Standard  Deviation 
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TABLE  IV- 5 


MEAN  17-41  SECOND  RMS  NOISE  AMPLITUDES  IN  m/x 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

9.  73 

3.  61 

8.  80 

3.  67 

10.  01 

3.  32 

CHTO 

8.  46 

3.  21 

9.  45 

4.  14 

9.  96 

3.  98 

MAIO 

7.  70 

2.  96 

7.  80 

2.  73 

8.  07 

% 

2.  91 

CTAO 

9.  41 

3.  29 

9.  90 

3.  02 

8.  78 

3.  18 

ZOBO 

7.  15 

2.  53 

7.  87 

2.  43 

7.  40 

2.21 

KAAO 

9.  40 

3.  78 

12.  53 

7.  59 

9.  61 

4.  96 

MAJO 

6.42 

1.  56 

7.  26 

1.  90 

7.  75 

2.  05 

S.  D.  = Standard  Deviation 


1 


IV-40 


1 


TABLE  IV- 6 
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MEAN  40-64  SECOND  RMS  NOISE  AMPLITUDES  IN  m/i. 


Vertical 

North 

East 

Old  11UU 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

3.  12 

1.08 

3.  23 

1. 08 

3.  95 

1.24 

CHTO 

2.  73 

0.  62 

3.  40 

0.  92 

3.  39 

1. 22 

MAIO 

3.  08 

2.  22 

3.  57 

1.  50 

4.  03 

% 

2.  32 

CTAO 

1.95 

0.  52 

3.  65 

1.  04 

3.66 

1. 48 

ZOBO 

2.  65 

1.15 

4.  24 

1.  94 

3.  56 

1. 09 

KAAO 

4.  03 

1.93 

9.46 

9.  28 

6.  06 

4.  66 

MAJO 

3.  92 

1.  57 

6.  00 

1.  70 

4.  74 

2.  05 

S.  D.  = Standard  Deviation 
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Figures  IV- 18a  through  IV- 23c  show  monthly  noise  level  trends 
which  are  derived  from  the  Julian  day  plots.  ANMO  10-25,  17-41,  and  40-64 
second  RMS  noise  trends  are  shown  in  Figures  IV-18a,  IV-18b,  and  IV-18c, 
respectively.  The  noise  level  for  ANMO  shows  a clear  pattern  as  it  drops 
from  January  through  June  and  rises  again  through  November.  A significant 
rise  in  the  noise  level  in  July  is  possibly  explained  by  thunderstorm  activity 
(Personal  Communication,  John  Hoffman,  1978).  While  trends  are  evident 
in  the  10-25  second  passband,  they  become  less  pronounced  in  the  17-41  sec- 
ond passband.  In  the  40-64  second  passband  no  noise  trend  can  be  discerned. 

Station  CTAO,  represented  in  Figures  IV-21a,  IV-21b,  and 
IV-21  c,  shows  a similar  phenomenon  between  passbands,  though  some  small 
trend  is  discernable  in  the  horizontal  components  in  the  40-64  second  pass- 
band.  This  trend  points  to  a noise  peak  in  September  and  October.  Other 
passbands  show  a clear  peak  in  July,  which  may  be  due  to  winter  storm 
activity  in  the  southern  hemisphere. 

MAIO,  in  Figures  IV-20a,  and  IV-20b,  representing  the  10-25 
and  17-41  second  passbands,  follows  a trend  similar  to  that  of  ANMO.  How- 
ever, the  40-64  second  passband  shows  a trend  which  does  not  appear  to  re- 
late to  the  traditional  seasons. 

KAAO  shows  noise  values  decreasing  from  July  through  Sep- 
tember in  Figures  IV-23a,  IV-23b,  and  IV-23c.  Values  rise  after  November 
in  the  1C-25  and  17-41  second  passbands,  but  continue  to  decline  in  the  40-64 
second  passband.  Unlike  stations  ANMO  and  CTAO,  KAAO  noise  trends  ap- 
pear most  pronounced  in  the  40-64  second  passband,  at  least  for  the  horizon- 
tal components. 

CHTO,  represented  in  Figures  TV-19a,  IV-19b,  andIV-19c, 
shows  a December  peak  in  the  10-25  and  17-41  second  trends.  This  peak  is 
absent  in  the  40-64  second  passband. 
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FIGURE  IV- 18c 

ANMO  40-64  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 1 9a 
CHTO  10-25  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 19b 

CHTO  17-41  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 19c 

CHTO  40-64  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 20a 

MAIO  10-25  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 20b 

MAIO  17-41  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 20c 

MAIO  40-64  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 21b 

CTAO  17-41  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 21c 

CTAO  40-64  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 22a 


ZOBO  10-25  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 22b 

ZOBO  17-41  SECOND  RMS  NOISE  TRENDS 
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ZOBO  40-64  SECOND  RMS  NOISE  TRENDS 
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FIGURE  IV- 23a 

KAAO  10-25  SECOND  RMS  NOISE  TRENDS 


FIGURE  IV- 23c 

KAAO  40-64  SECOND  RMS  NOISE  TRENDS 
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Figures  IV-24  through  IV-30  show  the  average  RMS  amplitude 
spectra  for  the  noise  at  each  station  under  evaluation.  These  spectra  are  un- 
corrected for  instrument  response.  They  were  determined  by  averaging  the 
individual  spectra  computed  from  the  noise  samples  and  converting  them  to 
RMS  amplitudes  using  Parseval's  formula: 


1 


where 


RMS 


Af 


a * \f‘  £ K>|2*C(£.) 


= the  elemental  frequency  interval  (Af  = 0.  C01953  Hz), 


A(f.f  a 


C(f.)  = 


the  discrete  Fourier  transform  spectral  density  esti- 
mate at  frequency  f.  , 

the  instrument  response  correction  at  frequency  f.  , 


a = the  initial  frequency  index,  and 

b = the  final  frequency  index. 

Since  the  RMS  amplitude  at  each  discrete  frequency  was  desired,  in  this  case 
a=b.  Also,  since  no  instrument  response  corrections  were  made,  C(f^)=l 
for  all  frequencies. 

The  left-hand  side  of  each  figure  shows  the  average  RMS  noise 
amplitude  spectra  for  each  component  of  motion.  These  spectra  were  com- 
puted for  the  period  range  of  11.6  to  85.  3 seconds  so  that  the  analyst  can  in- 
spect the  noise  levels  in  and  around  the  signal  periods.  The  right-hand  side 
of  each  figure  shows  the  log  mean  spectra  with  vertical  bars  representing 
plus  and  minus  one  standard  deviation.  These  are  presented  to  give  a quali- 
tative idea  of  the  day-to-day  variability  of  the  noise  from  period-to-period. 

The  following  comments  can  be  made  about  the  individual  fig- 
ures: 
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AVERAGE  RMS  AMPLITUDE  SPECTRA  ANMO  LONG- PERIOD  NOISE 
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AVERAGE  RMS  AMPLITUDE  SPECTRA  CHTO  LONG- PERIOD  NOISE 


SITE  50 


AVERAGE  RMS  AMPLITUDE  SPECTRA  CTAO  LONG- PERIOD  NOISE 


• Figures  IV-24  and  IV-25  - ANMO  and  CHTO  noise;  The  three 
components  show  nearly  identical  spectral  content  and  level  of 
the  local  noise  field.  At  each  station  the  noise  field  has  its 
greatest  day-to-day  variability  between  14  and  25  seconds. 

• Figure  IV-26  - MAIO  noise;  The  relative  height  of  the  17- 
second  microseismic  peaks  are  small  compared  to  the  rest  of 
the  MAIO  spectra.  At  periods  below  30  seconds,  the  three 
components  show  nearly  identical  spectral  content  and  level  of 
the  local  noise  field.  The  horizontal  components  show  a sec- 
ondary peak  at  37  seconds  nearly  as  large  as  the  microseismic 
peak.  The  noise  field  has  its  greatest  variability  between  14 
and  23  seconds  period  for  all  components  and  at  37  seconds 
period  for  the  horizontal  components. 

• Figure  IV-27  - CTAO  noise;  At  periods  below  30-seconds  the 

three  components  show  very  similar  spectral  content  and  level 

% 

of  the  local  noise  field.  The  vertical  component  17- second 
microseismic  peak  seems  to  be  truncated.  At  periods  above 
30  seconds  the  horizontal  components  show  a strong  rise  in 
noise  level  while  vertical  component  noise  level  increases  only 
slightly. 

• Figure  IV- 28  - ZOBO  noise;  At  periods  below  30  seconds  the 
three  components  show  nearly  identical  spectral  content  and 
level  of  the  local  noise  field.  At  periods  above  30  seconds, 
the  components  separate  as  horizontal  component  noise  levels 
rise  slightly  above  vertical  component  noise,  which  remains 
at  low  levels. 

• Figure  IV-29  - KAAO  noise;  At  periods  below  25  seconds  the 
three  components  show  similar  spectral  content  and  level  of 
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• I 

the  noise  field.  At  periods  above  25  seconds  the  horizontal 
components  show  increasing  noise  levels,  while  vertical  com- 
ponent noise  shows  little  increase  below  periods  of  about  70 
seconds. 

• Figure  IV-30  - MAJO  noise:  The  three  components  show  simi- 

lar spectral  content  and  noise  level  over  their  entire  length. 

Horizontal  component  noise  levels  rise  slightly  above  the  ver- 
tical component  noise  level  at  periods  above  50  seconds.  The 
17-second  microseismic  peaks  appears  rather  stunted  in  ampli- 
tude when  compared  with  the  rest  of  the  MAJO  spectra. 

In  general,  spectral  standard  deviations  seem  to  increase  with  increasing 
spectral  means.  Also,  no  evidence  was  found  from  these  spectra  supporting 
the  hypothesis  (Strauss  and  Weltman,  1977)  which  links  oceanic  wave  energy 
with  spectral  content  and  noise  level. 

Peak  20,  25,  and  30  second  noise  amplitudes  were  measured 
on  all  noise  samples  as  part  of  the  noise  analysis.  Each  value  measured 
represents  the  absolute  value  of  the  largest  zero-to-peak  excursion  of  a wave- 
form of  given  period  plus-or-minus  two  seconds.  Tables  IV-7,  IV-8,  and 
IV-9  list  the  means  and  standard  deviations  of  peak  and  log  peak  20,  25,  and 
30  second  measured  amplitudes. 
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Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

25.  59 

10.  90 

23.  52 

10.  44 

26.  04 

10.  04 

CHTO 

22.  17 

8.  16 

25.  42 

12.  19 

23.  95 

10.  60 

MAIO 

18.  13 

8.  55 

15.79 

7.  43 

15.  98 

9.  17 

CTAO 

26.27 

9.  97 

25.  74 

9 .01 

22.13 

9.  09 

ZOBO 

18.  00 

7.  35 

17.  99 

7.  51 

17.  57 

8.  04 

KAAO 

21.  33 

10.  34 

21.86 

12. 17 

20.  11 

10.45 

MAJO 

14.  71 

4.  59 

15.  64 

4.  62 

19.  35 

7.  87 

TABLE  IV- 7b 

MEAN  LOG10  PEAK  20  SECOND  NOISE  AMPLITUDES  IN  m/x 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

• 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

1.  37 

0.  17 

1.  33 

0.  19 

1.  38 

0.  17 

CHTO 

1.  31 

0.  18 

1.  36 

0.  19 

1.  33 

0.  19 

MAIO 

1.21 

0.  20 

1. 15 

0.  20 

1.  14 

0.  22 

CTAO 

1.  39 

0.  16 

1.  38 

0.  16 

1.  31 

0. 18 

ZOBO 

1.  22 

0.  17 

1.21 

0.  19 

1.  20 

0.  18 

KAAO 

1. 27 

0.  22 

1.27 

0.  23 

1. 24 

0.  22 

MAJO 

1.14 

0.  13 

1.  17 

0.  13 

1.  25 

0.  19 

S.  D.  = Standard  Deviation 


TABLE  IV- 8a 

MEAN  PEAK  25  SECOND  NOISE  AMPLITUDES  IN  m/x 


Station 

Vertical 

North 

East 

Mean 

I 

Mean 

S.  ] 

D. 

Mean 

S. 

D. 

ANMO 

22.  85 

10. 

05 

20. 

91 

10. 

24 

23.  78 

9. 

11 

CHTO 

19.  06 

7. 

54 

22. 

28 

8. 

42 

24.  96 

9. 

52 

MAIO 

20.  06 

7. 

45 

20. 

25 

7. 

76 

20.80 

8. 

39 

CTAO 

24.  65 

58. 

98 

21. 

82 

6. 

79 

24.  26 

35. 

08 

ZOBO 

17.  53 

6. 

85 

21. 

46 

8. 

30 

19.  84 

5. 

67 

KAAO 

23. 17 

9. 

43 

29. 

16 

13. 

55 

24.  32 

9. 

28 

MAJO 

17.44 

5. 

14 

19. 

50 

6. 

08 

21. 21 

5. 

50 

Standard  Deviation 
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TABLE  IV- 8b 

MEAN  LOGio  PEAK  25  SECOND  NOISE  AMPLITUDES  IN  m/i 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

1.  33 

0.  16 

1. 28 

0.  17 

1.  35 

0.  16 

CHTO 

1. 25 

0.  17 

1.  32 

0.  16 

1.  37 

0.  16 

MAIO 

1. 28 

0. 15 

1. 28 

0.  16 

1.  29 
% 

0.  16 

CTAO 

1. 26 

0.  22 

1. 32 

0. 13 

1.  30 

0.  20 

ZOBO 

1.  22 

0.  14 

1.  30 

0.  15 

1. 28 

0.  11 

KAAO 

1.  33 

0.  17 

1.42 

0.  20 

1.  36 

0.  16 

MAJO 

1. 22 

0.  12 

1. 27 

0.  12 

1.  31 

0.  11 

S.  D.  = Standard  Deviation 
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TABLE  IV- 9a 

MEAN  PEAK  30  SECOND  NOISE  AMPLITUDES  IN  m^t 


Vertical 

North 

East 

Station 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

12.  30 

6.  26 

12.  04 

6.  94 

14.  66 

6.  85 

CHTO 

9.  87 

4.  92 

14.  05 

8.  28 

13.  00 

5.  77 

MAIO 

13.  78 

6.  86 

17.  18 

8.  38 

18.  45 

8.  31 

CTAO 

11.  97 

18.  77 

12.  53 

6.  23 

13.  30 

7.  47 

ZOBO 

1 1 . 21 

5.14 

14.42 

5.  75 

12.  99 

4.  97 

KAAO 

16.  05 

8.  36 

29.  64 

24. 16 

21. 61 

23.  30 

MAJO 

10.  52 

4.  36 

15.99 

6. 17 

14.  22 

5.  05 

Standard  Deviation 
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TABLE  IV- 9b 

MEAN  LOG1q  PEAK  30  SECOND  NOISE  AMPLITUDES  IN  m/x 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

1. 03 

0.  22 

1. 02 

0.  24 

1.  11 

0.  21 

CHTO 

0.  94 

0.  22 

1. 09 

0.  21 

1.  07 

0.  19 

MAIO 

1. 09 

0.  20 

1. 09 

0.  21 

1. 07 
% 

0. 19 

CTAO 

0.  93 

0.  32 

1. 05 

0.  21 

1. 06 

0.  24 

ZOBO 

1. 00 

0.  22 

1.13 

0.  16 

1. 07 

0.  22 

KAAO 

1.  15 

0.  21 

1.  39 

0.  25 

1.  24 

0.  24 

MAJO 

0.  96 

0.  28 

1.  17 

0.  16 

1.10 

0.  25 

S.  D.  = Standard  Deviation 


i 


IV-76 


4 


4 


SECTION  V 

SRO  DETECTION  CAPABILITY 


A.  DISCUSSION 

The  Seismic  Research  Observatory  detection  capability  statis- 
tics for  the  five  new  stations  are  presented  in  this  section.  For  any  given 
event,  only  one  of  seven  conditions,  corresponding  to  seven  detection  codes, 
can  exist.  The  conditions  are  the  following: 

• Event  is  detected. 

• Event  is  not  detected. 

• Event  is  mixed. 

• No  data  are  recorded  for  the  time  period. 


Equipment  is  malfunctioning. 

Event  is  detected  on  microfiche  only. 

Event  epicenter  is  located  > 103°  from  station. 


Short-  Period 
Only 


A mixed  event  is  one  that  is  partially  or  completely  masked  by 
a second  signal.  This  happens  when  two  events  arrive  at  a station  at  essen- 
tially the  same  time  or,  when  a larger  signal  arrives  before  the  event  of 
interest,  burying  the  event  in  its  coda.  The  cause  of  no  data  being  recorded 
is  the  shutting  down  of  the  station.  Malfunctions  refer  to  the  partial  failure 
of  the  system,  from  the  sensor  unit  to  the  reception  of  data  at  the  Seismic 
Data  Analysis  Center,  which  causes  degradation  of  the  seismic  data.  A micro- 
fiche detection  only  refers  to  short-period  detection  capability,  and  will  be 
discussed  in  Subsection  V-B.  Since  epicentral  distances  greater  than  103° 
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constitute  the  P wave  shadow  zone,events  located  farther  than  103°  from  the 
station  of  interest  were  deleted  from  the  short-period  data  base  for  that  sta- 
tion. 

In  its  simplest  form,  a station's  detection  capability  would  be 
determined  by  whether  the  analyst  saw  either  the  event  of  interest  or  seismic 
noise.  However,  the  true  detection  capability  depends  greatly  on  how  mixed 
events,  malfunctions,  periods  of  no  recorded  data,  etc.  , are  treated.  For 
this  reason,  the  Seismic  Research  Observatory  detection  capability  estimates 
are  calculated  in  two  ways. 

The  first  of  these  is  labeled  the  'ideal  detection  capability1. 
When  calculating  this  estimate,  mixed  events,  events  for  which  no  data  was 
recorded,  and  events  containing  malfunctions  were  dropped  from  the  data 
base.  The  value  of  this  ideal  estimate  is  that  it  shows  the  detection  capability 
improvement  possible  if  the  reliability  of  the  instrumentation  could  be  im- 
proved and  if  methods  of  separating  mixed  events  could  be  found. 

The  second  estimate  is  labeled  'actual  detection  capability', 
and  it  considers  mixed  events,  events  for  which  no  data  was  recorded,  and 
events  containing  malfunctions  as  non-detections.  This  approach  gives  a 
real  world  detection  capability  estimate. 

The  quantity  used  to  represent  detection  capability  is  the  50 
percent  detection  threshold,  denoted  by  The  50  percent  detection 

threshold  is  the  bodywave  magnitude  for  which  the  probability  of  detection  is 
0.  5.  It  is  computed  by  fitting  the  Gaussian  probability  function  to  the  detec- 
tion statistics  by  a maximum  likelihood  method  (Ringdal,  1974). 

P SHORT- PERIOD  DETECTION  CAPABILITY  ESTIMATES 

In  the  previous  evaluations,  the  Seismic  Research  Observatory 
short- period  detection  capability  estimate  was  computed  from  data  recorded 


L 
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by  the  short-period  automatic  detector.  Therefore,  the  estimate  is  depen- 
dent upon  the  performance  of  the  detector.  In  an  effort  to  circumvent  this 
dependence,  the  daily  recordings  from  each  of  the  new  stations  were  obtain- 
ed on  microfiche  for  the  period  1 September  1977  to  30  November  1977,  and 
detection  statistics  were  generated  using  the  microfiche  as  data. 

The  criteria  for  determining  whether  an  event  was  detected  on 
the  microfiche  were  as  follows: 

• The  waveform  is  at  least  3.  5 dB  above  the  surrounding  noise 
waveform. 

• The  waveform  begins  within  ± 20  seconds  of  the  predicted  ar- 
rival time. 

The  requirement  that  the  maximum  amplitude  of  the  waveform 
be  at  least  3.  5 dB  above  the  noise  seems  low,  but  this  is  an  advantage  gained 
by  using  microfiche  as  compared  to  Calcomp  plots.  Rather  than  inferring 
the  noise  level  from  a small  segment  of  data;  noise  trends  could  be  seen,  and 
identification  of  an  event  could  be  made  with  more  confidence. 

Using  the  microfiche  also  enabled  the  analyst  to  pick  emergent 
waveforms;  however,  this  was  done  only  when  the  observed  start  time  satis- 
fied the  second  requirement.  It  should  be  noted  that  although  these  emergent 
start  times  were  chosen  at  the  most  obvious  break  from  the  noise  level,  they 
are  not  necessarily  the  actual  P wave  arrival  times,  but  could  be  as  much  as 
several  seconds  past  it. 

The  second  detection  criterion  is  not  absolute.  It  is  possible 
that  errors  in  computation  of  origin  time  and  location,  and  choice  of  start 
time  may  all  combine  to  place  the  first  point  of  detection  outside  the  ± 20 
second  gate.  Also,  the  expected  arrival  times  were  calculated  assuming  a 
normal  (33  km)  depth  of  focus.  In  cases  where  a P wave  was  observed  out- 
side the  gate,  and  no  other  events  could  be  found  in  the  National  Earthquake 
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Information  Service  (NEIS)  event  list,  a detection  was  declared.  These  cases 
were,  however,  very  rare. 
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In  order  to  provide  detection  codes  corresponding  to  the  condi- 
tions mentioned  at  the  beginning  of  this  subsection,  the  microfiche  detection 
statistics  were  then  compared  to  the  automatic  detector  logs  (see  Section  HI). 

This  generated  two  sets  of  detection  statistics,  one  for  the  analyst  detections 
and  one  for  the  automatic  detections.  Each  of  these  two  sets  were  further 
divided  in  the  manner  mentioned  in  the  first  part  of  this  section,  or,  into 
'ideal'  and  'actual'  detection  statistics.  The  four  resultant  detection  capabil- 
ity estimates,  calculated  from  these  four  sets,  are  shown  for  four  of  the  new 
stations  in  Figures  V-l  to  V-16.  Station  ZOBO  (Zongo,  Bolivia)  is  not  shown 
because  only  two  events  were  detected.  The  reason  for  this  was  too  large  an 
epicentral  distance  to  this  station  for  most  Eurasian  events. 

The  figures  labeled  'Ideal  Detection  Capability'  are  composed 
of  non- mixed  events  detected  both  by  the  automatic  detector  and  on  the  micro- 
fiche and  those  events  not  detected  by  either  one.  All  other  events  were  drop- 
ped from  the  detection  statistics,  i.  e.  , mixed  events,  events  for  which  no 
data  were  recorded,  events  containing  malfunctions,  and  events  detected  only 
on  the  microfiche.  This  is  the  ideal  detection  capability  as  'seen'  by  the 
automatic  detector. 

The  figures  labeled  'Actual  Detection  Capability'  are  composed 

o 

of  all  the  events  in  the  data  base  whose  epicenters  are  less  than  103  from 
the  station  of  interest.  However,  all  mixed  events,  events  for  which  no  data 
were  recorded,  events  containing  malfunctions,  and  events  detected  only  on 
thr  microfiche  were  considered  to  be  non- detections.  This  is  the  real-world 
c ection  capability  as  'seen'  by  the  automatic  detector. 

The  figures  labeled  'Ideal  Detection  Capability  - Detector  Dis- 
abled' are  identical  to  the  'Ideal  Detection  Capability'  except  that  the  former  accept 
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events  detected  only  on  the  microfiche  as  detections.  This  is  the  ideal  detec- 
tion capability  that  would  be  computed  by  an  analyst  viewing  continuous  data. 

The  figures  labeled  'Predicted  Detection  Capability  - Detector 
Disabled'  are  identical  to  the  'Actual  Detection  Capability'  except  that  the  former 
accept  events  detected  only  on  the  microfiche  as  detections.  This  is  a predic- 
tion of  the  actual  detection  capability  that  would  have  been  calculated  by  an 
analyst  viewing  continuous  data.  The  difference  between  the  'Actual  Detection 
Capability'  and  the  'Predicted  Detection  Capability  - Detector  Disabled'  is  the 
loss  in  detection  capability  caused  by  use  of  the  automatic  detector. 

The  upper  portion  of  each  figure  is  a histogram  giving  the  de- 
tection statistics  as  a function  of  bodywave  magnitude  for  that  particular  case. 
The  lower  portion  shows  the  percentage  of  events  detected  at  each  bodywave 
magnitude  (represented  by  asterisks),  the  fitted  maximum  likelihood  curve 
(represented  by  a solid  line),  and  the  90  percent  confidence  limits  for  this 
curve  (represented  by  dashed  lines).  The  values  for  'MB50'  and  'MB90' 
shown  on  the  figures  are  the  50  and  90  percent  detection  thresholds  as  picked 
from  the  maximum  likelihood  curve.  The  value  shown  for  'SIGMA'  i6  the 
standard  deviation  of  the  Gaussian  probability  function  obtained  by  the  maxi- 
mum likelihood  method. 


The  titles  for  each  figure  define  the  station  under  evaluation 
and  the  type  of  estimation,  as  defined  previously,  from  which  the  detection 
statistics  were  derived.  The  short-period  detection  capability  estimates  de- 
rived from  Figures  V-l  to  V-16  are  summarized  in  Table  V-l. 

The  MAJO  and  CHTO  ideal  and  actual  detection  capability  esti- 
mates differ  the  most  among  the  four  stations,  while  the  CTAO  estimates  differ 


the  least.  The  automatic  detector,  on  the  average,  cost  the  stations  0.  12  m^ 

units  of  detection  capability.  However,  the  detection  loss  for  CHTO  is  0.  20 

m,  units  and  for  MAJO  0.  15  m units,  two  to  three  times  the  loss  observed 
b b 

for  both  CTAO  and  KAAO. 
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TABLE  V-l 


SRO  SP  DETECTION  CAPABILITY 


50  Percent  Detection  Threshold 

Calculation 

Type 

CHTO 

CTAO 

KAAO 

MAJO 

mb 

S.  D. 

mb 

S.  D. 

mb 

S.  D. 

mb 

S.  D. 

Ideal 

4.  60 

0.  06 

5.43 

0.  10 

4.  26 

0.  07 

4.  91 

0.  09 

Actual 

4.  95 

0.  08 

5.  51 

0.  11 

4.  45 

0.  07 

5.  29 

0.  12 

Ideal 

(Detector 

Disabled) 

4.  53 

0.  06 

5.  38 

0.  09 

4.  24 

0.  07 

4.  82 

0.  08 

Predicted 

(Detector 

Disabled) 

4.  75 

' 

0.  06 

5.  45 

0.  10 

4.  37 

0.  08 

5.  14 

0.  10 

S.  D.  = Standard  Deviation 


A 


c. 


LONG- PERIOD  DETECTION  CAPABILITY  ESTIMATES 
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Estimates  of  long-period  Seismic  Research  Observatory  detec- 
tion capability  are  presented  for  each  of  the  five  stations  evaluated  in  this  re- 
port. The  criteria  for  determining  whether  a detection  has  been  achieved  for 
a given  event  are; 

• The  presence  of  dispersion  in  the  signal  gate. 

• The  presence  of  a peak  in  the  dispersed  wave  train  3 dB  or 
more  above  any  peak  outside  the  dispersed  wave  train  and  in- 
side a time  gate  starting  600  seconds  before  the  predicted 
Love  wave  arrival  time  and  ending  600  seconds  after  the  esti- 
mated Rayleigh  wave  end  time. 

• The  occurrence  of  the  signal  onset  within  ±180  seconds  of  the 
predicted  signal  onset  time. 

• Detection  of  the  event  on  at  least  two  components. 

% 

Occasionally,  an  event  is  considered  to  be  detected  when  not  all  criteria  are 
satisfied.  Signal  peaks  are  occasionally  less  than  3 dB  above  the  noise  peaks 
and  still  recognized  as  signals  from  their  dispersion  characteristics.  Also, 
at  a given  station  it  is  sometimes  possible  to  find  specific  features  of  a seis- 
mic waveform  from  a given  region,  which  enables  the  analyst  to  detect  the 
event  even  though  not  all  detection  criteria  have  been  satisfied.  An  interest- 
ing example  of  this  is  discussed  in  the  preliminary  Seismic  Research  Observa- 
tory report  (Strauss,  1976).  The  last  criterion  is  intended  to  reduce  the  num- 
ber of  false  alarms,  i.  e.  , declaring  an  event  detected  when  it  is  not. 

The  problem  of  mixed  events  is  often  difficult  for  the  analyst 
to  resolve  and  may  be  a source  of  false  alarms.  (The  term  'false  alarm'  in 
this  case  means  declaring  a detection  when  in  fact  the  observed  signal  is  from 
an  event  other  than  that  under  analysis.  ) When  a signal  is  observed  in  the 
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time  gate  of  the  event  under  analysis,  the  analyst  first  checks  the  waveforms 


on  the  three  components  of  motion  to  see  that  their  inter-relationships  are 
correct.  If  still  in  doubt,  the  analyst  checks  available  event  lists  to  see 
whether  any  other  reported  event  could  have  arrived  in  the  signal  gate.  In 
general,  the  analyst  declares  a detection  if  a dispersed  signal  is  observed 
having  the  correct  inter-relationships  between  the  Love  and  Rayleigh  waves 
and  if  no  other  event  has  been  reported  which  could  be  mistaken  for  the  event 
under  analysis. 

The  long- period  detection  statistics  and  derived  maximum 
likelihood  curves  are  presented  in  Figures  V-17  to  V-26.  The  detection 
capability  estimates  were  calculated  in  the  manner  described  in  the  first  por- 
tion of  this  section,  i.  e.  , an  'ideal  detection  capability'  and  an  'actual  detec- 
tion capability'. 

Figures  V-17  to  V-26  are  displayed  in  the  same  manner  as 
Figures  V-l  to  V-16  for  the  short-period  detection  capability  estimates,  ex- 
cept that  microfiche  were  not  used  for  the  long-period  detection  statistics. 
Since  no  detector  is  involved,  there  are  only  two  long-period  detection  capabil- 
ity estimates  for  each  station  rather  than  four,  as  for  short-period.  Each  of 
the  figures  represents  one  long-period  detection  capability  estimate,  where 
the  upper  portion  shows  a histogram  of  the  detection  statistics  and  the  lower 
portion  shows  the  maximum  likelihood  curve  fitted  to  these  statistics. 


The  long-period  detection  capability  estimates  derived  from 
these  figures  are  summarized  in  Table  V-2.  There  is  considerable  difference 
between  the  ideal  and  actual  estimates  for  all  the  stations  except  CTAO.  One 
possible  reason  for  the  over-all  higher  difference,  as  compared  to  the  short- 
period  detection  capabilities,  is  the  greater  possibility  for  long-period  events 
tc  -e  mixed.  ZOBO  and  MAJO  have  especially  large  differences;  ZOBO  be- 
cause of  the  greater  distance  and  thus  greater  chance  for  a mixed  event,  and 
MAJO  because  of  the  large  amount  of  system  failure  (resulting  in  no  data  be- 
ing recorded). 
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50  Percent  Detection  Threshold 


Calculation 

Type 

CHTO 

CTAO 

KAAO 

MAJO 

ZOBO 

mb 

S.  D. 

S.  D. 

“b 

S.  D. 

mb 

S.  D. 

mb 

S.  D. 

Ideal 

4.  56 

0.  06 

5.  05 

0.  08 

4.  07 

0.  09 

4.  36 

0.  09 

5.  26 

■ 

Actual 

4.  90 

0.08 

5.  20 

0.  08 

4.  53 

0.  07 

4.  97 

0.  08 

5.  91 

0.  26 

SECTION  VI 

EARTHQUAKE- PRESUMED  EXPLOSION  DISCRIMINATION 


A.  DISCUSSION 

This  section  considers  the  question  of  discriminating  between 
earthquakes  and  presumed  nuclear  explosions  using  long-period  data.  The 
discrimination  method  used  is  the  surface  wave  magnitude  (M^)  versus  body- 
wave  magnitude  (m^)  plot.  This  plot  of  Mg  versus  mfa  is  expected  to  function 
as  an  earthquake -presumed  nuclear  explosion  discriminant,  since  for  a given 
m^,  an  explosion  produces  much  weaker  Rayleigh  waves  than  does  an  earth- 
quake. Weichert  and  Basham  (1973)  assign  the  Mg-  m^  criterion  a 30  per- 
cent probability  of  correctly  identifying  an  explosion,  accompanied  by  a false 
alarm  probability  of  about  0.  1 percent. 


B.  COMPUTATION  OF  SURFACE  WAVE  MAGNITUDES 


The  equations  used  in  Mg  computation  are  adapted  by  Strauss 
(Personal  Communication,  1978)  from  data  presented  by  Nutley  and  Kim 
(1975).  Mr.  Strauss  and  the  authors  feel  that  the  Mg-  m^  curve  is  best  re- 
solved into  three  segments,  based  on  epicentral  distance: 
for  A S 10°  Mg  = Log1()A/T  + 0.  89  Log1()A  + 1.  34 

10°<A<25°  M = Log.-A/T  + 1.  07  Log  A + 1.16 
8 10  10 

A > 25°  Mg  = Logj  qA/T  + 1. 66  Log1()A  + 0.  34 

where 


A = epicentral  distance  in  degrees 
A = peak-to-peak  amplitude  in  millimicrons 
T = period  in  seconds  of  the  measured  amplitude. 


M values  were  calculated  automatically  from  identified  detected  events 
s 

stored  on  event  tapes  (Section  II).  The  Mg  computation  routine  chooses  an 
edit  gate  based  on  expected  surface  wave  arrival  times  and  durations  stored 
in  each  event  header.  The  times  of  zero  crossings  within  the  edit  gate  are 
calculated,  and  the  maximum  absolute  amplitudes  in  millimicrons  between 
adjacent  zero  crossings  are  found.  The  associated  waveform  periods  are 
calculated  as  two  times  the  time  between  these  adjacent  zero  crossings. 
Epicentral  distances  and  sample  rates  are  read  from  the  event  header.  Quan- 
tization factors  and  instrument  response  corrections  are  built  into  the  pro- 
gram. (Instrument  response  corrections  are  derived  from  values  supplied 
by  Albuquerque  Seismological  Laboratory  via  AFTAC.  ) Mg  values  are  cal- 
culated and  printed  out.  The  largest  Mg  value  at  each  specified  period  of 
interest  is  displayed  separately.  The  analyst  may  now  check  these  largest 

M values  against  the  plot  of  the  data.  If  he  does  not  like  an  M value  at  a 
s 

particular  period  (for  example,  the  Mg  value  appears  to  be  associated  with 
a noise  pulse  in  the  signal  gate)  he  can  select  another  value  from  the  Mg  list 
generated  by  the  program. 


\ ' 


C.  DISCRIMINATION  RESULTS 

Table  VI- 1 lists  those  events  of  the  data  base  which  may  be 
termed  presumed  nuclear  explosions.  They  were  selected  on  the  basis  of 
their  epicentral  locations,  bodywave  magnitudes,  and  shallow  depths.  Thus, 
the  events  from  Flinn  and  Engdahl  Region  41  (Flinn  and  Engdahl,  1965)  were 
selected  because  their  epicenters  are  at  the  Nevada  Test  Site  and  their  body- 
wave  magnitudes  are  much  larger  than  known  earthquakes  in  the  area.  The 
events  from  Region  329  were  selected  because  their  epicenters  are  at  the 
e tern  Kazakh  test  site  and  their  bodywave  magnitudes  are  larger  than  most 
earthquakes  from  that  area. 
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♦Detection  Status  Codes 

1 - Event  detected  according  to  detection  criteria  of  Section  V 

2 - Event  not  detected 

3 - Mixed  event 

4 - System  failure,  no  data  recorded 

5 - Malfunction  (glitch,  spike,  or  the  like)  in  signal  gate. 


The  numbers  in  the  columns  under  the  station  headings  in 
Table  VI- 1 refer  to  the  detection  status  of  the  events  as  recorded  at  each  sta- 
tion. The  number  1 indicates  that  the  event  was  detected,  2 indicates  that  it 
was  not  detected,  3 indicates  that  it  was  mixed  with  another  event,  4 indicates 
no  data  were  recorded  for  that  event,  and  5 indicates  that  a system  malfunc- 
tion obscured  the  data.  These  codes  indicate  that  six  of  the  fourteen  events 
cannot  be  used  in  the  discrimination  study,  since  they  were  not  detected  at 
any  of  the  stations. 

Table  VI- 2 presents  M values  from  the  various  stations  which 

s 

detected  presumed  nuclear  explosions.  The  average  difference  between  M 

s 

and  m^  is  1.  6 units.  Hudson  and  Douglas  (1975)  estimated  that  for  a given 

m , an  explosion  should  yield  surface  wave  magnitudes  about  0.  5 M units 
b s 

lower. 

Figures  VI-1  through  VI- 3 show  vertical  component  25- second 
surface  wave  magnitudes  for  stations  CHTO,  CTAO,  ZOBO,  KAAO,  and 
MAJO.  Symbols  used  are:  o - earthquake;  • - presumed  explosion.  All  pre- 
sumed nuclear  explosions  are  also  labeled  by  their  event  numbers.  The 
straight  line  in  each  plot  represents  the  Mg-  relationship  for  that  data  set 
computed  using  all  earthquake  data  points. 


Referring  to  Figures  VI-1  through  VI-3  the  following  comments 


can  be  made; 


Event  1805  shows  no  separation  from  either  the  CHTO,  CTAO, 

or  MAJO  earthquake  Mg  population  (Figure  VI- 1 and  Figure 

VI- 3)  although  the  M values  for  that  event  lie  on  the  lower 

s 

bound  of  each  concerned  station's  earthquake  population.  Ac- 
cording to  the  Mg-  discriminant  this  event  may  have  been 
an  explosion  or  an  earthquake. 
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TABLE  VI- 2 


CHTO  AND  CTAO  V 


ZOBO  AND  KAAO  VERTICAL  COMPONENT  25-SECOND  M - m DATA 


VERTICAL  COMPONENT  25 -SECOND  M - m,  DATA 


• Event  1874' s Mg  value  is  well  separated  from  CHTO's  and 
ZOBO's  earthquake  Mg  populations  (Figure  VI-1).  This  dis- 
criminant would  classify  it  as  an  explosion. 

• Events  1630  and  1640  yielded  Mg  values  which  are  well  separated 
from  the  KAAO  earthquake  population  (Figure  VI- 2)  but  only 
marginally  separated  from  the  earthquake  population  of  MAJO 
(Figure  VI- 3).  These  events  also  appear  to  have  been  explo- 
sions. 

• Events  1772  and  1773  yielded  Mg  values  which  separate  mar- 
ginally from  earthquake  M populations  at  MAJO  (Figure  VI- 3). 

s 

The  M value  recorded  at  station  KAAO  for  event  1772  also 

8 

shows  marginal  separation  from  earthquakes  recorded  at  that 
station  (Figure  VI-2).  These  events  are  also  probably  explo- 
sions but  more  uncertainty  exists  here  than  in  the  classification 
of  the  three  previously  discussed  events. 

• Event  1693's  Mg  value  lies  on  the  periphery  of  the  KAAO  earth- 
quake population  (Figure  VI- 2).  Classification  of  this  event  is 
inconclusive. 

• Event  1781' s Mg  value  lies  well  within  MAJO's  earthquake  popu- 
lation (Figure  VI-3).  This  event  classifies  as  an  earthquake. 

In  summary,  while  no  station  is  infallible  when  considered 
alone,  CHTO,  CTAO,  ZOBO,  KAAO,  and  MAJO  show  potential  as  contributors 
to  an  Mg-  m^  discrimination  network. 

D.  Mg-  mb  RELATIONSHIPS 

Table  VI- 3 presents  the  results  of  fitting  a straight  line  to  the 
Mg-  m^  points.  The  Mg-  m^  relationships  are  presented  only  for  the  case  of 
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TABLE  VI- 3 


SRO  STATIONS  Ms-  mb  RELATIONSHIPS  FOR  VERTICAL 
COMPONENT  25-SECOND  ENERGY 


Station 


CHTO 

CTAO 

ZOBO 

KAAO 

MAJO 


where: 


0.  08 
0.07 
0.  04 
0.  08 


Center  of  Mass 

“b 

M 

s 

n 

4.  87 

4. 12 

83 

4.  99 

4.  33 

50 

4.  89 

4.  78 

19 

4.  81 

4.  25 

93 

4.  79 

4.  11 

66 

-4.  92 
-4.  31 
-2.66 
-3.  71 
■21.  30 


= amb  + b 

= variance  normal  to  the  Mg-m^  estimate, 
and 

= number  of  Mg-m^  points  used  in  compu- 
tation. 
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all  data  points  measured  at  each  station.  Note  the  relatively  steep  slope  and 
large  negative  intercept  of  the  MAJO  Mfl-  rr^  line.  No  explanation  has  been 
found. 

The  reader  is  referred  to  Figures  VI- 1,  VI- 2,  and  VI- 3 for 
pictorial  representations  of  both  the  M#-  relationships  presented  here 
and  the  data  points  with  which  these  relationships  were  established. 
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SECTION  vn 
CONCLUSIONS 


This  section  summarizes  the  results  of  this  continued  evalua- 
tion of  the  Seismic  Research  Observatories.  The  major  results  are: 

A.  DATA  QUALITY 

• In  general,  the  data  quality  was  good.  At  all  stations  except 
ZOBO,  data  were  degraded  by  malfunctions  less  than  4%  of 
the  time.  Most  ZOBO  malfunctions  were  caused  by  a defective 
field  tape.  Other  malfunctions  were  mostly  glitches  in  the  sig- 
nal gate. 

• Estimates  of  station  reliability  were  based  on  the  frequency  of 
occurrence  of  malfunctions  and  on  estimates  of  station  down- 
time. ZOBO  and  MAJO  both  showed  low  reliabilities  of  ap- 
proximately 0.  7 due  to  relatively  large  amounts  of  station 
down-time.  The  other  stations  had  reliability  figures  above 

0.  8. 

• The  frequency  of  mixed  event  occurrence  for  long-period  wave- 
forms ranged  from  a low  of  0. 19. for  CHTO  and  MAJO  to  a high 
of  0.  46  for  station  ZOBO.  The  frequency  of  mixed  event  oc- 
currence is  high  for  ZOBO  Eurasian  events  because  of  the  long 
distance  between  that  station  and  Eurasia  as  compared  to  other 
stations. 


VII- 1 


B.  SHORT- PERIOD  AUTOMATIC  DETECTOR 


• During  the  time  period  that  the  detector  was  studied,  detection 
probability,  Pp,  ranged  from  50%  to  80%  while  the  false  alarm 
rate  (FAR)  varied  between  0.  2 and  2.  1 false  alarms  per  hour. 

P seemed  rather  independent  of  FAR;  this  may  reflect  FAR 
instability  inherent  in  straight  STA/LTA-type  detectors. 

• Detector  parameter  setting  and  freezing  practices  have  im- 
proved from  1976  to  1977  as  evidenced  by  a more  stable  rela- 
tionship between  percent  time  recorded  and  number  of  detec- 
tions. 

C.  NOISE  ANALYSIS 

• In  order  of  increasing  magnitude, the  mean  short-period  RMS 
noise  values  in  the  0.  5-4.  0 Hz  passband  (uncorrected  for  in- 
strument response)  are:  1.02  m n at  ZOBO,  1.  67  m/z  at  CHTO, 
1. 94  mfi  at  KAAO,  3.  53  m/x  at  MAJO,  and  5.  55  m fi  at  CTAO. 
These  values  show  a strong  correlation  with  station-nearest 
coastline  separation,  indicating  that  a majority  of  short-period 
noise  is  derived  from  ocean  wave  energy  injected  at  the  coast- 
line. 

• Long-period  RMS  noise  level  and  trends  were  investigated  in 
the  passband s:  10-25,  17-41,  and  40-64  seconds.  A year's 
worth  of  noise  data  were  available  for  three  stations,  ANMO, 
MAIO,  and  CTAO.  Noise  levels  for  these  stations  seemed  to 
change  seasonally  in  the  10-25  and  17-41  second  passbands. 
However,  trends  in  the  40-64  second  passband  noise  levels 
are  not  easily  explained. 
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• In  general,  no  correlation  was  suggested  by  the  data  between 
long-period  noise  level  or  spectral  content  and  station- 
coastline  separation. 

• A comparison  between  VLPE  and  SRO/ASRO  instrument  re- 
sponse corrected  spectra  for  the  same  sites  show  significant 
improvements  in  noise  levels  after  the  installation  of  borehole 
seismometers.  ASRO  instruments  showed  strong  correlation 
between  noise  amplitude  and  noise  variance  in  contrast  to  their 
VLPE  counterparts. 

DETECTION  CAPABILITY 

• Actual  50%  m^  short-period  Eurasian  detection  capabilities 
are:  4.  95  at  CHTO,  5.  43  at  CTAO,  4.  26  at  KAAO,  and  4.  91 
at  MAJO. 

• Only  two  short-period  events  were  detected  by  ZOBO,  proving 
this  station  a poor  detector  of  Eurasian  events. 

• Actual  50%  m^  long-period  Eurasian  detection  capabilities  are: 
4.  56  at  CHTO,  5.  05  at  CTAO,  4.  07  at  KAAO,  4.  36  at  MAJO, 
and  5.  26  at  ZOBO. 

• Mixed  events  had  very  little  effect  on  short-period  detection 
capabilities.  Their  effect  was  significant  on  long-period  detec- 
tion capability,  especially  at  ZOBO. 

• The  automatic  short-period  detector  as  compared  to  analyst 
microfiche  detection,  showed  an  average  loss  of  0.  1 2 units 
in  detection  capability. 

• Station  down-time  cost  MAJO  about  0.  5 m^  units  in  detection 
capability.  Down-time  for  ZOBO  also  significantly  affected 
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its  detection  capability.  Actual  detection  capabilities  for  both 
stations  should  improve  since  all  apparent  hardware  problems 
have  been  solved. 


E.  DISCRIMINATION 

• Stations  CHTO,  CTAO,  ZOBO,  KAAO,  and  MAJO  classified 

four  out  of  fourteen  presumed  nuclear  explosions  with  some 
confidence  on  the  basis  of  an  Mg-m^  criterion.  Some  indica- 
tion of  classification  was  noted  for  all  but  one  of  the  events 
which  were  detected,  suggesting  that  these  stations  would 
function  best  as  part  of  a discrimination  network. 


The  following 

Mg-  m^  relationships  were  established. 

CHTO  : 

M 

s 

= 1 . 86  - 4.  92 

CTAO  : 

M 

= 1.  73  m - 4.  31 

s 

b 

ZOBO  : 

M 

s 

= 1.  52  - 2.  66 

KAAO  : 

M 

s 

= 1. 66  - 3.  71,  and 

MAJO  : 

M 

s 

= 5.  33  - 21.  30 

The  anomalous  relationship  for  MAJO  could  not  be  explained. 

Future  work  on  the  evaluation  of  the  Seismic  Research  Obser- 
vatory stations  should  be  directed  toward  the  following  points: 

• Evaluate  all  stations  as  they  become  operational. 

• Expand  the  short-period  and  long-period  noise  data  bases  so 
that  at  least  one  year  of  data  is  available  at  each  station  for 
the  study  of  noise  level  trends. 

• Data  bases  should  be  individually  tailored  for  at  least  those 
stations  which  are  far  from  Eurasia.  In  this  respect,  ZOBO 
should  be  re-evaluated  in  such  a way  that  its  detection  capabil- 
ity is  measured  more  relevantly. 
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APPENDIX  A 


THE  DATA  BASE 


The  247  events  used  for  this  evaluation  of  the  Seismic  Research 
Observatory  stations  are  listed  on  the  following  pages.  Short-period  detection 
capability  data  bases  included  only  those  events  with  station- epic  enter  dis- 
tances of  less  than  103°. 

In  the  data  base  list,  the  column  headed  'EVNO'  gives  the  iden- 
tifying number  assigned  to  each  event.  The  columns  headed  'DATE'  and 
'TIME'  give  the  date  and  origin  time  of  each  event.  The  epicentral  coordin- 
ates are  listed  under  'LAT.  ' and  'LONG.  ' for  latitude  in  degrees  north  and 
longitude  in  degrees  east.  The  bodywave  magnitude  is  listed  under  the  head- 
ing 'MB'.  The  column  headed  'Q'  gives  the  NORSAR  quality  rating  for  the 
parameters  of  each  event  where 

1 = good  to  excellent  quality 

2 = fair  to  good  quality 

3 = poor  to  fair  quality. 

The  seismic  source  area  is  listed  in  the  column  headed 
'LOCATION'.  Finally,  the  geographic  region  as  defined  by  Flinn  and  Engdahl 
(Flinn  and  Engdahl,  1965)  for  each  event  is  listed  under  the  heading  'SUBREG'. 
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APPENDIX  B 

INSTRUMENT- RESPONSE  CORRECTED  NOISE  CHARACTERISTICS 


In  the  course  of  the  evaluation  of  the  Seismic  Research  Obser- 
vatory stations,  instrument- response  corrected  noise  field  characteristics 
were  studied  at  each  station  under  evaluation.  In  this  appendix  to  the  con- 
tinuation of  the  Seismic  Research  Observatories  evaluation,  the  extension  to 
this  study  is  presented. 

To  correct  the  data  for  the  effects  of  the  instrument  response, 
the  data  were  first  Fourier  transformed  into  the  frequency  domain.  Noise 
amplitude  spectra  were  then  computed  and  smoothed  to  64  frequencies.  The 
smoothing  was  performed  to  simplify  the  task  of  making  the  instrument  re- 
sponse corrections.  In  the  course  of  this  smoothing,  the  spectral  values 
were  scaled  so  that  Parseval's  formula  would  hold.  This  formula,  which  per- 
mits the  computation  of  RMS  noise  from  frequency  domain  data,  is  expressed 


where 


A(£i>r 


[Af  53  I A(f.)|  2 *C(f.)  | (frequency 

. i=a  1 1 J estima 


aency-domain 

stimate) 


= the  elemental  frequency  interval, 

= the  discrete  Fourier  transform  spectral  density 
estimate  at  frequency  f , 

= the  instrument  response  correction  at  frequency  f , 
= the  initial  frequency  index,  and 
= the  final  frequency  index. 


$ 


Figures  B-la  through  B-6c  show  the  instrument- response  cor- 
rected RMS  noise  values  plotted  against  Julian  day  for  the  passbands  of  10-25, 
17-41,  and  40-64  seconds.  Figures  B-7a  through  B-12c  show  the  RMS  noise 
level  trends  computed  from  these  values.  Since  these  figures  differ  from 
those  of  Section  IV,  where  the  instrument  response  was  not  corrected  for, 
only  in  level,  no  further  discussion  will  be  presented. 

Tables  B-l  through  B-3  present  means  and  standard  deviations 
of  instrument  response  corrected  noise  in  the  10-25,  17-41,  and  40-64  second 
passbands.  Tables  B-4  and  B-5  present  means  and  standard  deviations  of  in- 
strument response  corrected  peak  and  log  peak  20-  and  30- second  noise  am- 
plitudes. Instrument  response  corrected  mean  peak  and  mean  log  peak  ampli- 
tudes were  derived  from  their  non-instrument  response  corrected  analogues 
through  the  implementation  of  the  following  statistical  formulae: 

Given  x 

x 

Log  x 

where  x'  is  a peak  noise  amplitude  and  c is  an  instrument  response  correc- 
tion factor.  Instrument  response  corrected  and  non- instrument  response  cor- 
rected 25 -second  values  are  equivalent  since  the  SRO  long-period  response  is 
normalized  at  that  period.  The  reader  is  again  referred  to  Section  IV  for 
discussion. 

Figures  B-lb  through  B-l 9 present  the  mean  instrument  - 
rf  ponse-corrected  RMS  amplitude  spectra  for  each  of  the  seven  stations 
evaluated  here.  The  right  half  of  each  figure  displays  Log^  ^ RMS  amplitude 
spectra.  The  associated  vertical  bars  represent  plus-or-minus  one  standard 
deviation.  These  stations,  with  the  sole  exception  of  MAJO  (Figure  B-l 9) 
display  a microseismic  peak  at  a period  of  about  17  seconds  followed  by  a 


= x'/c 


1 . 1 
— x’  ; cr  = — cr 
c x c x' 


= Log  x'  - Log  c ; 


= cr , 
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ANMO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND  RMS  NOISE 


SITE  30 

Vertical 


ANMO  INSTRUMENT  RESPONSE  CORRECTED  40-64  SECOND  RMS  NOISE 


rtical 


CHTO  INSTRUMENT  RESPONSE  CORRECTED  10-25  SECOND  RMS  NOISE 


MAIO  INSTRUMENT  RESPONSE  CORRECTED  10-25  SECOND  RMS  NOISE 
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FIGURE  B-3b 

MAIO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND  RMS  NOISE 


SITE  36 

Vertical 


MAIO  INSTRUMENT  RESPONSE  CORRECTED  40-64  SECOND  RMS  NOISE 


CTAO  INSTRUMENT  RESPONSE  CORRECTED  10-25  SECOND  RMS  NOISE 


SITE  50 

Vertical 


CTAO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND  RMS  NOISE 


B-14 


4 


1 


SITE  52 

Vertical 


KAAO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND  RMS  NOISE 
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FIGURE  B-7a 

ANMO  INSTRUMENT  RESPONSE  CORRECTED  10-25  SECOND 

RMS  NOISE  TRENDS 


FIGURE  B- 7b 

ANMO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND 

RMS  NOISE  TRENDS 
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Vertical 


FIGURE  B-7c 

ANMO  INSTRUMENT  RESPONSE  CORRECTED  40-64  SECOND 

RMS  NOISE  TRENDS 
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FIGURE  B-  8a 

CHTO  INSTRUMENT  RESPONSE  CORRECTED  10-25  SECOND 

RMS  NOISE  TRENDS 


FIGURE  B-8b 

CHTO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND 

RMS  NOISE  TRENDS 
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FIGURE  B-  9b 

MAIO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND 

RMS  NOISE  TRENDS 
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SITE  36 

Vertical 


FIGURE  B-9c 

MAIO  INSTRUMENT  RESPONSE  CORRECTED  40-64  SECOND 

RMS  NOISE  TRENDS 
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FIGURE  B-lOb 


CTAO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND 

RMS  NOISE  TRENDS 
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FIGURE  B-lOc 

CTAO  INSTRUMENT  RESPONSE  CORRECTED  40-64  SECOND 

RMS  NOISE  TRENDS 
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FIGURE  B-llb 


ZOBO  INSTRUMENT  RESPONSE  CORRECTED  17-41  SECOND 

RMS  NOISE  TRENDS 
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FIGURE  B-llc 

ZOBO  INSTRUMENT  RESPONSE  CORRECTED  40-64  SECOND 

RMS  NOISE  TRENDS 
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FIGURE  B-12a 

KAAO  INSTRUMENT  RESPONSE  CORRECTED  10-25  SECOND 

RMS  NOISE  TRENDS 


FIGURE  B-  12c 

KAAO  INSTRUMENT  RESPONSE  CORRECTED  40-64  SECOND 

RMS  NOISE  TRENDS 
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TABLE  B-2 

INSTRUMENT  RESPONSE  CORRECTED  MEAN  17-41  SECOND 
RMS  NOISE  AMPLITUDES  IN  m/a 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

12.  67 

4.  63 

11. 40 

4.  83 

12.  07 

4.  37 

CHTO 

11.12 

4.  39 

12.  14 

5.  64 

12.93 

5.  36 

MAIO 

9.  47 

3.  70 

9.45 

3.  33 

9.  77 

3.  70 

CTAO 

12.  44 

4.  44 

13.  00 

4.  11 

11. 42 

4.  36 

ZOBO 

9.  07 

3.  25 

9.  91 

3.  33 

9.  29 

3.  00 

KAAO 

11.  50 

4.  72 

14.  96 

8.  50 

11. 65 

5.  81 

MAJO 

7.  85 

1. 95 

8.  76 

2.  37 

9.  48 

2.  68 

S.  D.  = Standard  Deviation 


TABLE  B-3 

INSTRUMENT  RESPONSE  CORRECTED  MEAN  40-64  SECOND 
RMS  NOISE  AMPLITUDES  IN  m/x 


Station 

Vertical 

No  rth 

East 

Mean 

ESI 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

5.  64 

2.  10 

5.  74 

1. 98 

7.  10 

2.  23 

CHTO 

4.  96 

1.  22 

6.  02 

1.  55 

6.  06 

2.  36 

MAIO 

5.46 

4.  11 

6.  27 

2.  71 

4.  39 

CTAO 

3.46 

0.  96 

6.  65 

2.  81 

6.  63 

1.  99 

ZOBO 

4.  69 

2.  04 

7.  58 

3.  77 

6.  36 

1.  87 

KAAO 

7.  14 

3.47 

11. 49 

7.  85 

11.13 

8.  88 

MAJO 

7.  06 

8.  69 

3.  79 

8.  62 

3.  92 

S.  D.  = Standard  Deviation 


B-41 


4 


TABLE  B-4a 


INSTRUMENT  RESPONSE  CORRECTED  MEAN  PEAK  20  SECOND 
RMS  NOISE  AMPLITUDES  IN  m/x 


Station 

Vertical 

North 

East 

Mean 

S. 

D. 

Mean 

S. 

D. 

S. 

D. 

ANMO 

31. 

13 

13. 

26 

28. 

61 

12. 

70 

31.  68 

12. 

21 

CHTO 

26. 

97 

10. 

48 

30. 

92 

14. 

83 

29.  14 

12. 

89 

MAIO 

22. 

06 

10. 

40 

19. 

21 

9. 

04 

19.  44 

11. 

15 

CTAO 

31. 

96 

12. 

13 

31. 

31 

10. 

95 

26.  92 

11. 

06 

ZOBO 

21. 

90 

8. 

94 

21. 

89 

9. 

14 

21. 37 

9. 

78 

KAAO 

25. 

95 

12. 

58 

26. 

59 

14. 

80 

24.46 

12. 

71 

MAJO 

17. 

89 

5. 

59 

19. 

03 

5. 

62 

23.  54 

9. 

57 

S.  D.  = Standard  Deviation 
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TABLE  B-4b 

INSTRUMENT  RESPONSE  CORRECTED  MEAN  LOG  PEAK 
20  SECOND  NOISE  AMPLITUDES  IN  m/x 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

1. 46 

0.  17 

1.42 

m 

1. 47 

0.  17 

CHTO 

1. 40 

0.  18 

1.45 

0.  19 

1. 42 

0.  19 

MAIO 

1.  30 

0.  20 

1. 24 

0.  20 

1.  23 

0.  22 

CTAO 

1.48 

0.  16 

1.47 

0.  16 

1. 40 

0.  18 

ZOBO 

1. 31 

0.  17 

1.  30 

0.  19 

1.  29 

0.  18 

KAAO 

1.  36 

0.  22 

1.  36 

0.  23 

1.  33 

0.  22 

MAJO 

1.  23 

m 

1.  26 

0. 13 

1.  34 

0.  19 

S.  D.  = Standard  Deviation 


B-43 


TABLE  B- 5a 


■i 

t- 


INSTRUMENT  RESPONSE  CORRECTED  MEAN  PEAK  30  SECOND 
NOISE  AMPLITUDES  IN  m/z 


Station 

Vertical 

North 

East 

Mean 

S,  D. 

Mean 

S.  D. 

Mean 

S. 

D. 

ANMO 

12. 

50 

6.  36 

12. 

24 

7.  05 

14.  90 

6. 

96 

CHTO 

10. 

03 

5.  00 

14. 

28 

8.  41 

13.  21 

5. 

86 

MAIO 

14. 

00 

6.  97 

17. 

46 

8.  52 

18.  75 

8. 

45 

CTAO 

12. 

16 

19.  08 

12. 

73 

6.  33 

13.  52 

7. 

59 

ZOBO 

11. 

39 

5.  22 

14. 

65 

5.  84 

13.  20 

5. 

05 

KAAO 

16. 

31 

8.  50 

30. 

12 

24.  55 

21.  96 

23. 

68 

MAJO 

10. 

69 

4.  43 

16. 

25 

6.  27 

14.  45 

5. 

13 

B -44 


* 


TABLE  B-5b 

INSTRUMENT  RESPONSE  CORRECTED  MEAN  LOG 

30  SECOND  NOISE  AMPLITUDES  IN  m/x  10 


Station 

Vertical 

North 

East 

Mean 

S.  D. 

Mean 

S.  D. 

Mean 

S.  D. 

ANMO 

1. 04 

0.  22 

1.  03 

0.  24 

1.  12 

0.  21 

CHTO 

0.  95 

0.  22 

1.  10 

0.  21 

1. 08 

0. 19 

MAIO 

1.  10 

0.  20 

1.18 

0.  24 

1.  23 

0.  20 

CTAO 

0.  94 

0.  32 

1.  06 

0.  21 

1.  07 

0.  24 

ZOBO 

1.  01 

0.  22 

1.  14 

0.  16 

1. 08 

0.  22 

KAAO 

1.16 

0.  21 

1. 40 

0.  25 

1. 25 

0.  24 

MAJO 

0.  97 

0.  28 

1.18 

0.  16 

1.  11 

0.  25 

S.  D.  = Standard  Deviation 
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minimum  in  the  period  range  25-28  seconds.  Noise  amplitudes  beyond  this 
minimum  appear  to  rise  almost  monotonically.  This  spectral  structure  is 
also  evident  at  NWAO  and  SNZO,  and  to  some  extent  at  GUMO  and  TATO 
(Strauss  and  Weltman,  1977).  Similar  spectral  structures  have  been  publish- 
ed for  ambient  earth  motion  measurements  at  Queens  Creek,  Arizona;  LASA; 
Las  Cruces,  New  Mexico;  and  Ogdensburg,  New  Jersey  (Fix,  1972).  Spectral 
analyses  at  Murphy  Dome,  Alaska  (Von  Seggern,  1970);  and  of  the  TFSO  ar- 
ray (Chiburis,  1967)  also  indicated  strong  microseismic  noise  peaks  at  17 
seconds.  The  lack  of  a spectral  peak  at  station  MAJO  may  be  due  to  the  fact 
that  only  twenty  noise  samples  are  involved  in  the  average. 

Reprinted  for  comparison  purposes  in  Figures  B-20  through 
B-24  are  noise  spectra  generated  for  the  VLPE  stations  located  at  or  near 
SRO/ASRO  stations,  see  Table  B-6  (Prahl,  1974).  The  Chiang  Mai  and 
Albuquerque  VLPE  stations  were  replaced  with  SRO  borehole  seismometers. 
As  was  expected,  noise  amplitudes  dropped  significantly.  We  furthermore 
observe  that,  in  contrast  to  the  VLPE  station  noise,  the  spectral  amplitude 
variances  of  SRO/ASRO  station  noise  are  strongly  correlated  with  the  noise 
spectral  amplitudes.  At  CTAO,the  variances  are  less  strongly  correlated 
than  other  SRO/ASRO  stations  (Figure  B-16).  Also,  the  spectral  amplitudes 
and  amplitude  variances  are  larger  at  VLPE  CTA  than  at  ASRO  CTAO.  A 
comparison  of  SRO  ZOBO  spectra  with  those  of  VLPE  ZLP  reveal  markedly 
distinct  characteristics.  ZOBO's  three  components  appear  practically  iden- 
tical below  30  seconds  and  remain  similar  at  higher  periods.  ZLP's  three 
components  show  much  less  similarity  to  one  another.  However,  these  two 
stations  do  not  share  the  same  location  as  do  the  other  station  pairs,  which 
may  account  for  the  spectral  differences  observed.  On  the  other  hand,  SRO 
station  MAJO  and  VLPE  station  MAT  do  share  locations  and  these  two  stations 
also  show  distinctly  different  spectra.  While  the  three  component  interrela- 
tionships seem  normal  for  both  sites,  the  MAJO  spectra  appear  damped  be- 
tween periods  of  12  and  20  seconds.  It  should  be  remembered,  though,  that 
the  MAJO  spectra  are  preliminary  and  based  on  only  20  noise  samples. 
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ENT  STATION  CMC 


B-56 


4 


ENT  STATION  MAT 


TABLE  B-6 

VERY  LONG  PERIOD  EXPERIMENT  STATIONS  AND  LOCATIONS 


Station  Name 

Station 

Designator 

Numeric 

Designator 

Latitude 

Longitude 

Charters  Towers, 
Australia 

CTA 

1 

20. 09S 

146.  26E 

Chiang  Mai, 

Thailand 

CHG 

2 

18.  79N 

98.  98E 

Albuquerque, 

New  Mexico 

ALQ 

9 

34.  94N 

106.  46  W 

La  Paz, 

Bolivia 

ZLP 

10 

16. 50S 

68. 13W 

Matsushiro, 

Japan 

MAT 

11 

36.  54N 

138.  21 E 

